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RICHARD MOECKEL AND RICHARD MONTGOMERY 



Abstract. We carry out a sequence of coordinate changes for the planar 
three-body problem which successively eliminate the translation and rotation 
symmetries, regularize all three double collision singularities and blow-up the 
triple collision. Parametrizing the configurations by the three relative position 
vectors maintains the symmetry among the masses and simplifies the regu- 
larization of binary collisions. Using size and shape coordinates facilitates the 
reduction by rotations and the blow-up of triple collision while emphasizing the 
role of the shape sphere. By using homogeneous coordinates to describe Hamil- 
tonian systems whose configurations spaces are spheres or projective spaces, 
we are able to take a modern, global approach to these familiar problems. We 
also show how to obtain the reduced and regularized differential equations in 
several convenient local coordinates systems. 



1. Introduction and History 

The three-body problem of Newton has symnietries and singularities. The re- 
duction process eliminates symmetries thereby reducing the number of degrees of 
freedom. The Levi-Civita regularization eliminates binary collision singularities by 
a non-invertible coordinate change together with a time reparametrization. The 
McGehee blow-up eliminates the triple collision singularity by an ingenious polar 
coordinate change and another time reparametrization. Each process has been ap- 
plied individually and in various combinations to the three-body problem, many 
times. 

In this paper we apply all three processes globally and systematically, with no 
one body singled out in the various transformations. The end result is a complete 
flow on a five-dimensional manifold with boundary. We focus attention on the 
geometry of the various spaces and maps appearing along the way. At the heart 
of this paper is a beautiful degree-4 octahedral covering map of the shape sphere, 
branched over the binary collision points (see figure |4]). This map first appears in 
the work of Lemaitre [TH [13] ■ One of our goals is to give a modern, geometrical 
approach to this regularizing map. 

The reduction procedure for the three body problem dates back to Lagrange 
who found elegant differential equations for 10 translation and rotation invariant 
variables, including the the squares of the lengths of the three sides of the triangle 
formed by the bodies. These equations are valid for the three-body problem in any 
dimension. The variables of Lagrange also have the advantage of maintaining the 
symmetry among the masses. On the other hand, for the planar problem they are 
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subject to 3 nonlinear constraints in addition to the energy and angular momentum 
integrals. Moreover, we do not know a way to regularize the binary collision singu- 
larities in Lagrange's equations. For a modern introduction to Lagrange's equations 
see [lllig. 

Jacobi eliminates the translation symmetry by the familiar device of fixing the 
center of mass at the origin and introducing Jacobi coordinates [9j . The elimination 
of rotations is achieved by introducing some angular variable (or variables in the 
spatial case) to describe the overall rotation of the triangle together with some com- 
plementary, rotation-invariant variables. This method, which is the basis for much 
of the later work on the three-body problem, has some disadvantages. First, the 
Jacobi coordinates break the symmetry among the masses, making it much more 
difficult to regularize all three binary collisions at once. Second, for topological 
reasons, there is no way to choose an angular variable suitable for a global reduc- 
tion which includes the binary collision configurations. Namely, the map from the 
normalized configuration space to the shape sphere is a Hopf fibration, a nontrivial 
circle bundle. If we delete the binary collision points, the bundle becomes trivial 
but this deletion is not conducive to subsequent regularization. 

Murnaghan f24| derived a symmetrical Hamiltonian for the planar three-body 
problem in terms of the lengths of the sides and an angular variable representing 
the overall rotation of the triangle with respect to an inertial coordinate system. 
Then he obtains a reduced system by ignoring the angular variable. Van Kampen 
and Wintner carry out a similar reduction for the spatial three-body problem [29j . 
While these reductions avoid breaking the symmetry, they are still subject to the 
problem about the use of angular variables in a nontrivial bundle. In addition, using 
the side lengths as variables leads to differential equations which are not smooth at 
the collinear configurations (a problem seemingly avoided somehow by Lagrange). 

Lemaitre 112] introduced a symmetrical approach to reduction and regularization 
of binary collisions leading to the octahedral branched covering map of the sphere 
mentioned above. After using Euler angles to reduce by rotations, he introduces a 
size variable and two shape variables which can be viewed as spherical coordinates 
on the shape sphere which we use below. The regularization of binary collisions is 
done in the shape variables by means of the octahedral covering map. The use of 
Euler angles limits the validity of the reduction step of Lemaitre's work and the 
derivations are based on rather heavy trigonometric computations. But much of 
this paper can be viewed as a modern, global way to arrive at his covering map. 

In this endeavor we have the advantage of the modern theory of reduction of 
Hamiltonian systems with symmetry. Smale describes the reduction process for the 
three-body problem as the formation of a quotient manifold with a reduced Hamil- 
tonian flow [271 ■ Meyer [T5] and Marsden-Weinstein [TB| formalized the reduction 
procedure into what is now called "symplectic reduction theory" . Fixing the inte- 
grals of motion determines invariant manifolds in phase space. The quotient spaces 
of these invariant manifolds are the reduced phase spaces and the flows induced on 
them are again Hamiltonian with respect to an appropriate symplectic structure 
and a reduced Hamiltonian function. 

The regularization procedure goes back to Levi-Civita [141 who showed how 
to regularize binary collisions in perturbed planar Kepler problems by using the 
complex squaring map (a branched double covering of the complex plane). It is 
easy to adapt his method to regularize one of the binary collisions in the three-body 
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problem, but regularizing all three requires more ingenuity. Lemaitre's regularizing 
map behaves like the complex squaring map at each of the binary collision points 
on the shape sphere. Another approach to simultaneous regularization (without 
reduction) was introduced by Waldvogel [30 who used a quadratic mapping of the 
translation-reduced configuration space C^. We use a similar quadratic mapping 
applied to certain homogeneous shape variables below. Heggie fS' found an elegant, 
symmetrical way to regularize all of the binary collisions for the A^-body problem. 
In the planar case, his method is to apply separate Levi-Civita transformations to 
each of the difference vectors qi — qj. We apply this same idea below, but to the 
homogeneous shape variables, where it is found to induce Lemaitre's octahedral 
covering. 

Triple collision acts like an essential singularity in the three-body problem. 
McGehee in 1974 [T^ showed how an extension of spherical coordinates, together 
with a time reparametrization yields a fiow with no singularities at triple collision. 
This "McGehee blow-up" has the effect of replacing the triple collision point by 
a manifold called the collision manifold. Relative to the new parameterization, it 
takes forever to reach triple collision, whereas the Newtonian time to triple collision 
is finite. The fiow on the triple collision manifold governs the behavior of near-triple 
collision solutions. One aspect of the blow-up procedure is the use of separate size 
and shape coordinates to describe the configuration of the bodies. As shown below, 
such a splitting also facilitates the global reduction by rotations. 

Several authors have combined blow-up of triple collision with reduction and/or 
regularization of binary collision. Waldvogel reduced and regularized the fiow on 
the zero-angular-momentum triple collision manifold [31] . The first part of his 
paper combines Murnaghan's reduction procedure with some formulas of Lemaitre 
to obtain a reduced and regularized Hamiltonian for the zero-angular momentum 
three-body problem. Binary collisions are not regularized on the nonzero angular 
momentum levels. However, it is known that triple collisions can only occur when 
the angular momentum is zero. After restricting to the zero angular momentum 
manifold, Waldvogel blows up the triple collision to get reduced, regularized and 
blown-up differential equations. Simo and Susin used these coordinates in their 
study of the dynamics on the collision manifold |26]. These coordinates are very 
much in the spirit of this paper but do not achieve a full reduction, regularization 
and blow-up due to the restriction to zero angular momentum. 

The present paper draws on all these sources. We begin with some symplec- 
tic reduction theory. Turning to the three-body problem we eliminate translation 
symmetry by introducing the three difference vectors Qij = qi — qj as coordinates. 
Since these are linearly dependent, some effort is needed to justify the change of 
coordinates. Next we introduce a size variable and an associated spherical coordi- 
nates Xij. One novelty of our approach is that we use homogeneous coordinates 
to describe points on spheres. Instead of constraining the spherical coordinates to 
have a fixed norm, we only ask them to avoid the origin and then we find differential 
equations for them which are invariant under scaling. 

Once this point of view is adopted, it is relatively easy to carry out a global 
reduction by rotations. Using complex coordinates, the combined action of scaling 
and rotation is just scaling by a complex number. Quotienting by complex scaling 
we end up with a complex projective space, in fact with CP^. Of course as real 
manifolds CP"'^ ~ and this is our version of the shape sphere. We finally obtain a 
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global reduction of the planar three-body problem with a six- dimensional reduced 
phase space, the cotangent bundle of M+ x S^. 

Turning to regularization, we use simultaneous Levi-Civita transformations of 
the homogeneous variables Xij to regularize all three binary collisions. This regu- 
larizing map is applied to both the rotation-reduced and unreduced problems. In 
the reduced case we get a reduced and regularized system on the cotangent bundle 
of E"*" X which is related to the unregularized version by Lemaitre's map. 

Finally we show how McGehee's blow-up procedure can be applied to the various 
Hamiltonians we have found. 

2. Symplectic Reduction 

In this section we recall some results about the reduction of a Hamiltonian system 
with symmetry. In addition we show how to tell when two symmetric Hamiltonian 
systems lead to equivalent reduced systems. 

First we describe the basic symplectic reduction theory of Meyer [18] and Marsden- 
Weinstein [16' in the case of a system with symmetry. Suppose (M, a;) is a sym- 
plectic manifold and G is a Lie group which acts on M as a group of symplectic 
diffeomorphisms. Let J : Af — ^ g* be the momentum map, where q* is the dual 
of the Lie algebra of G. If we fix a momentum value /i S 0* and suppose that the 
action of G maps the level set J~^(/i) into itself. The quotient space 

p^ = J-\^i)/G 

is called the reduced phase space. 

If the group action is free and proper, then this space is a smooth manifold. 
There is an induced symplectic form lo^ on which is obtained as follows. First, 
for X € M, restrict ijj{x) to the tangent spaces TxJ~^{fJ,). The resulting two-form 
has a kernel which is precisely the tangent space to the group orbit through x. This 
implies that there is an induced two-form on the quotient vector space which is the 
tangent space to the quotient manifold. 

Now a H : M ^ R is a G-invariant Hamiltonian then the corresponding Hamil- 
tonian flow has J~^{^) as an invariant set and G-orbits map to G-orbits under 
the flow. Hence there is a well-defined quotient flow on J~^{ii)/G. There is also 
a reduced Hamiltonian iJ^ : — > K and the reduction theorem states that the 
quotient flow on (P^,a;^) is the Hamiltonian flow of the reduced Hamiltonian. 

Now suppose we have two such Hamiltonian systems with symmetry. For i — 1,2, 
there will be symplectic manifolds (Mi,uji), symmetry groups Gi and momentum 
maps Ji. If we fix momentum values fii we get reduced reduced phase spaces 
Pi = J^^{^i)/Gi with symplectic forms ui^-. Suppose Hi : Mi — >■ M are Gi-invariant 
Hamiltonians and let iJ^. : M be the corresponding reduced Hamiltonians. 

We want to give a concrete way to check that the two reduced Hamiltonian flows 
are equivalent. 

Suppose we have a smooth map F : J^^{^i) — > J2^{^J,2) which maps Gi-orbits 
into G2-orbits, i.e., F is equivariant. Then F induces a smooth map of quotient 
manifolds P : Pi — > P2. We will call F partially symplectic if it preserves the 
restrictions of the symplectic forms, i.e., F* {uj2\j-i(^^^^^) = wijj-i^^^^. It follows 

that F : (Pi,a;pj) — )■ (P2,a;^2) is symplectic. Hence P is a local diffeomorphism, 
even if P itself is locally neither injective nor surjective. Then the usual theory of 
symplectic maps applied to P gives: 
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Theorem 1. Suppose F : J^^{fii) — > J2^{fj,2) is a partially symplectic, equivariant 
map and that the restrictions of the Hamiltonians are related by Hi = H20F . Then 
F : Pi ^ P2 is a symplectic, local dijfeomorphism of the reduced phase spaces which 
takes orbits of the reduced Hamiltonian flow of Hf^-^ to those of H^^ . 

Definition 1. A partially symplectic, equivariant map G : ,12^(1^2) Ji^il^i) 
such that F o G — id(niodG2) and G o F — id(modGi) (so that these maps take 
group orbits into group orbits) will be called a pseudo-inverse for F . 

A partial inverse G for F induces a bona fide inverse G for F which exhibits an 
equivalence between the two reduced Hamiltonian flows. 

As a special case, suppose the two Hamiltonians are both defined on the same 
space and have the same symmetry group. If their restrictions to J^^(/i) agree 
then they will lead to the same reduced system. The identity map will provide the 
required partially symplectic map. We will call two such Hamiltonians equivalent. 
Equivalent Hamiltonians may produce different flows on J~^(/^) but the quotient 
flows will agree. 

The following theorems about the symplectic reduction of a cotangent bundle 
M — T*X will be used later. (See [5] Theorem 4.3.3 for a version of these theorems.) 
Suppose G acts freely on the configuration space X and that the G-action on M 
is the canonical lift of this base action. Suppose that the orbit space B for the G 
action on X is a manifold and the projection tt ; X — > _B a submersion. 

Theorem 2. Under the above assumptions, the reduced space Pq of T* X at fi ^ 
is isomorphic to T* B with its canonical symplectic structure ljb 

The theorem can be proved as a special case of Theorem [T] Because tt is onto, 
dux : TxX — > T^(^x)B is an onto linear map for each x & X. Consequently the dual 
map c?7r* : T*^^-^B — > T*X is injective. In the next paragraph we will show that the 

image of this dual is J^^{0)x'- 

(1) im{dTTl) = J-\Q)x := J-i(O) n T*X. 

It follows that we can invert d7r* on the fiber J^^{{))x C T*X. Define 

F : J-i(O) ^ T*B-F{x,p) = (7r(x), 

One verifies that F is a partially symplectic map relative to G acting on J~^(0), 
and the trivial group acting on T*B. A particularly easy way to see the partially 
symplectic nature of F is to introduce local bundle coordinates X D 7r~^(C/) = 
UxG. {X is covered by sets of this nature. ) In bundle coordinates Tr{x, g) — x and 
so T*X ^ T*UxGxQ*. We write elements o{T*X over U as (&, P; g, n),beU,P e 
T^U,g £ G,fj. £ Q*. In these coordinates J{b, P; g, ^) = fi, so that the general 
element of J~^{0)u can be written (6, Ph, g, 0) and F{b, Pt, g, 0) = (6, Pt). We have 
uJx = dx A dP + dg A d^ and, ivb — dx A dP, where we hope the meaning of these 
symbolic expressions is obvious. It follows immediately that F*lob = wjsf|,/-i(o) 
which is the claimed partially symplectic nature of F. Theorem [2] follows. 

We explain why eq ([!]) holds, and in the process gain some understanding of the 
momentum map. The group action is a map G x X —i' X which, when differentiated 
with respect to g S G at the identity yields the "infinitesimal action" a : q x X 
TX. For each frozen x, the map ■ Q ^ T^X is linear and, because G acts freely, 
injective. As we vary x, a forms a vector bundle map, part of an exact sequence of 
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vector bundle maps over X: 

^ Q-x X -^-^ TX n*TB 

where tt*TB = {{x,V);x & X,V & T^(x)B} is the puh-back of TB over B by 
the map tt : X ^ B. (Exactness of the sequence foUows by differentiating the 
statement that the fibers of tt are the G-orbits.) Duahzing we get 

0* X X T*X n*T*B 

The momentum map for the G-action on T*X is tti o a* where tti : g* x X — > g* 
is the projection onto the first factor. In other words 

J{x,p) = (J* p. 

It fohows from the exactness of the dual sequence that im{dTT*) — ker{a*) which 
is precisely equation ([T]). 

In order to identify the reduction of M = T*X at a non-zero value, /i ^ 0, 
we introduce a connection F for the bundle G ^ X ^ B. The curvature of the 
connection F is a g- valued two- form on B, which we may pull-back to T*B via 
the canonical projection tb : T* B ^ B. Then fi ■ Q is a, scalar- valued two-form on 
B. 

Theorem 3. Under the same assumptions as above on G, the reduced space of 
T* X at fi is isomorphic to T* B with the twisted symplectic structure ujb — r^/z • f2 

We only present the proof in the case G = S^, whose Lie algebra we identify 
with K in the usual way. Then a connection is a G-invariant one- form on T*X 
which satisfies the normalization property: 

j{x,r{x)) = 1. 

Its curvature is defined by 

dT = 

We define the momentum shift map: 

$^ : J-\0) ^ J-\p) = {x,p + pr{x)). 

which adds pT pointwise to each covector. The fiber-linearity of J shows that $^ 
does indeed map J'^iO) onto J'^{^). (The inverse of subtracts fxT. ) The 
map is G-equivariant since F is G-invariant. Thus $^ induces a G-equivariant 
diffeomorphim J~^(0)/G J^^{^i)/G. We have already identified J~-^(0)/G with 
T* B. However, $^ is not partially symplectic, so we cannot directly apply theorem 
1. To understand and quantify this failure, let G = PdQ denote the canonical 
one-form on T*X. Compute $* 9 = 9 + /ir^F. Taking the exterior derivative, 
using ujx — — c?9 we find that $^wx = — /iT^7r*f2. This equation implies that 
if we shift the canonical two- form on J^^(O) by subtracting iJLT'tn*^! then $p is a 
partially symplectic map between J~^(0) and J^^(/i). Theorem follows. 

3. Reduction by Translations 

To formulate the Newtonian planar three-body problem, it is convenient to use 
the complex plane where we identify (x, y) € with x + iy & C 

Let Qi S C,j = 1,2, 3 be the positions of the three-bodies and let q = ((?i, (72, 93) G 
C'^. We will adopt the Hamiltonian point of view where the conjugate momentum 
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variables pi are covectors rather than vectors. If we identify a covector (a, b) G 
M^* with a + ib G C then we have momentum variables G C* ~ C and p — 

bl,P2,P3) €C3*. 

The planar three-body problem is the Hamiltonian system on the phase space 
(C^ \ A) X C'''* with Hamiltonian 

H{q,p) = Koip)-Uiq) 



2mi 2m2 2m3 
mim2 m^mi m2mj, 



(2) ^ob) = 

I91-92I ks-gil 192-931 
where A = {g : = qj for some i ^ j}, the singular set. From now on, we will not 
explicitly mention that the singular set must be deleted from the domains of the 
various Hamiltonians we construct. 

The Newtonian potential is invariant under the group G — C acting by transla- 
tion on the position vectors and leaving the momenta fixed. The momentum map 
is given by 

Plot = Pi + P2 + J33 e C* . 
By fixing a value of this integral and passing to the quotient space one obtains 
a reduced Hamiltonian system. A simple and familiar way to accomplish this 
reduction is to assume ptot = and then fix the center of mass at the origin: 

miqi + 771292 + "73^3 = 0. 

However, we will now describe an alternative method for eliminating the trans- 
lation symmetry which will make it easier to regularize double collisions later on. 
This approach is a variation on the one used by Heggie in [8j. We will view it as 
an application of theorem [T] 

3.1. Relative coordinates. Introduce relative position variables (3i2, Qsi, Q23 G 
C and corresponding momentum variables P12, P^i, P2'i G C*. The relative coordi- 
nates are related to the positions variables qi by a linear map Q = Lq 

(3) L : ^ Qr2 = 9i - 92 Qzi =93-91 Q23 = 92 - 93- 

The dual map, which describes the pull-back of the relative momenta Pij to p space, 
is given by 

(4) L*:C3*->C3* Pl=Pl2-P31 P2^P23-Pl2 P3^P31-P23- 

We naturally have Qjt — —Qij and consequently Pji ~ —Pij so that eq Q can be 
written pi = YijPij, a form which extends to the A'^-body problem. 
The linear map L is neither one-to-one nor onto. Its kernel, 

ker L — {q : q — {c, c) for some c G = C}, 

is the subspace of translation symmetries in q-space. So its image 

W - imi = {Q : Q12 + g3i + Q23 = 0} 

is isomorphic to the quotient space of C"^ by translations. W is a complex subspace 
of C'^ with complex dimension two, or real dimension 4. We can define a map in 
the other direction, q — L\Q) 

fr,-, j-i "l2Qi2 - ^3(331 m-3(523 - ™lQl2 "^iQsi - "^2Q23 

(5) : qi = ,92 = ,93 = 

m m ni 
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maps C'^ onto 



W' = iniL^ = {q : niiqi + TO2'72 + "^s^a — 0} 



the zero-center of mass subspace and it is easy to check that the restrictions i|vv' 
and i^lw are inverses. 

For the dual map, we find that the kernel is generated by translations in P- 
momentum space 



ker i* — {P : P — {c, c, c) for some c £ C 
while the image is the zero-momentum subspace 

V = imL* ={p:pi+p2+Pn = 0}. 



The map 

it* 
(6) Lt* , 



-<3* 



TO2P1 



miP2 



m 



Psi 



TO1P3 - m^pi 



m 



P 



23 



m3P2 - TO2P3 

m 



maps 



onto 



V' = imL'* = {P : TO3P12 + ^2^31 -|- TO1P23 
and the restrictions L*|v' and i^*|v a-re inverses. 



0} 



Define a relative coordinate Hamiltonian on the {Q,P) phase space x C'^* by 

Hrel{Q,P)=K{P)-U{Q) 



(7) 



KiP) = KoiL*P) 



P 



311 



\P 



23 



P12P 



27711 2to2 

mim2 m^rrii m2m^ 
IQ12I ^ IQ31I ^ IQ23I ' 
so that 

(8) H{q,L*P)^Hrel{Lq,P). 

The kinetic energy can be written 



-t- 



|/'31 



231 



2ma 



(9) K(P) = Ip^BP 



V( 1 
nil 



s = 



7711 

(— + — )^ 

^ 77I3 7711 

— ^-I 

7773 





(— + — )^ 

\ m2 mg ' - 



where / denotes the 2x2 identity matrix. 



3.2. Equivalance to the translation-reduced three-body problem. We will 

now show that the reduction of the Hamiltonian system with Hamiltonian TJ^e; (Q, P^ 
by translations in momentum space is equivalent to the reduction of the three-body 
Hamiltonian iJ by translations in configuration space. 

Theorem 4. W x C"^* is invariant under the Hamiltonian flow of HreiiQ, P)- The 
restricted flow is invariant under translations in momentum space and it induces a 
quotient flow which is conjugate to the zero total momentum flow of the three-body 
problem reduced by translations. 

The proof will be an application of theorem[T] First we describe how the relevant 
symplectic structures look in complex coordinates. If Q G and P G C^* it is 
convenient to define a Hcrmitian variant of the natural evaluation pairing: 



(10) 



{P,Q)^Pl2Ql2 + P3lQ31 + 
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As a result, if Qjk — Xjk + i Ujk and Pjk = ajk + ihjt we get 

re(P, Q) = 012X12 + 612^12 + ■■ • 

im(P, Q) = ai22/i2 - 6122^12 + ■ • ■ 

Thus the real part of the complex pairing agrees with the usual real pairing and, 
as a bonus, the imaginary part is — /i where /i is the angular momentum. With 
this convention, the canonical one-forms on (g,p)-space and (Q,P)-space can be 
written 

e = re(p, dq) = re(pi dqi + p2 dq2 + ps dq^) 

e = re(P, dQ) = re(Pi2 dQi2 + P31 dQai + A3 ^^23)- 

Proof of theorem 2} For the three-body problem we have the phase space Mi = 
X C^* = {{q,p)} with the standard symplectic structure. The Hamiltonian 
H{q,p) is invariant under the action of the group Gi = C acting by c • {q,p) — 
{qi + c,q2 + c,q3 + c,pi,p2,P3),c G C. We fix the momentum level ptot = and 
obtain a quotient Hamiltonian flow. 

For the Hamiltonian Hrei the phase space is M2 = C"^ x C''* = {{Q, P)} with the 
standard symplectic structure. Hj-eiiQ, P) is invariant under the action of the group 
G2 = C* acting on by c • (Q, P) = (Q12, Q31 , Q23, P12 + c, P31 + c, P23 + c), c G C* . 
The momentum map is Qtot = Q12 + Q31 + <323 and we fix the momentum level 
Qtot = giving a second quotient Hamiltonian flow. 

To see that these two quotient flows are equivalent we apply Theorem [T] Define 

F{q,p) - {Lq, L^*p) G{Q, P) = (Ltg, L*P). 

Then, F : {ptot = 0} ^ {Qtot = 0} and G : {Qtot = 0} {ptot = 0}. Moreover 
GoF{q,p) = c - {q,p) where — c = ^(miqi +m2q2 + '^13(73) € C is the center of mass. 
Similarly F o G{q,p) = c • (Q,P) where -c = ;^(m3Pi2 + ™2f'3i + miP23) € C*. 
In other words G o F ~ id(mod Gi) and F o G — id(mod G2). 

It remains to verify that F and G are partially symplectic. Consider the canonical 
one-forms (12). From ([s]) and (|6|. We find, for example P*Pi2 ~ {jn2Pi — TniP2)/m 
and F*dQ 12 — dqi — dq2. After a bit of algebra we get 

P*e = 6 - re{ptotimidqi + m2dq2 + m^dq^i) / m) . 

Restricting to {ptot — 0} shows that F is partially symplectic. Similarly 

G*e = (d- re((TO3Pi2 + m2P3i + miP2z){dQi2 + dQ^i + dQ2z))lm 

which we restrict to {Qtot = 0} to see that G is also partially symplectic. We have 
shown that F and G are pseudo-inverses in the sense of definition [T] According to 
eq. ([s]) these pseudo-inverses intertwine H and H^ei ■ The hypothesis of theorem I 
have been verified, completing the proof. QED 

Hamilton's equations for the Hamiltonian HreiiQ, P) are simply 
Q = BP 

(13) i,_TT _ ,mim2Qi2 rn^rrnQsi TO2"13Q23s 

P -f^Q - -( 13 ' Z3 ' Z3 ) 

'12 '31 '23 

where r^j = \Qij\. (Note that here and in all of the differential equations below, 
partial derivatives like Uq are calculated by simply calculating the corresponding 
real partial derivatives and converting the resulting real vector or covector to com- 
plex notation; no complex differentiations are involved). Differential equations for 
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the three-body problem reduced by translations are obtained by restricting Q to 
W. Then Q remains in W under the flow. Moreover, covectors P, P' which are 
initially equivalent under translation remain so. 

Since the symmetry group C* acts only on the momenta P^ , the reduced phase 
space is the eight-dimensional space W x (C^*/C*) ~ W x iniL* W x V. This can 
be identified with the cotangent bundle T*W = W x W* as follows. Let P e C^*. 
Then Pjyy S W* and two covectors P, P' e C^* have the same restriction to W if 
they differ by an element of kerL*, i.e., if they are equivalent under the symmetry 
group. 

So far we have not really accomplished any "reduction" since there are still 12 
(Q, P) variables. Essentially, we have traded the constraint ptot = Pi +J'2 +P3 = 
and the translation symmetry in q for the constraint Qtot — Q12 + Q31 + Q23 — 
and translation symmetry in P. We will see below that the use of the Qij is advan- 
tageous for regularizing double collisions. A genuine reduction of dimension can be 
easily achieved by introducing a basis for W. Moreover, this can be accomplished 
in several ways as we will see in section 3.4 below. But one virtue of Q is that it 
avoids making a choice of parametrization and thereby preserves the symmetry of 
the problem under permutations of the masses. 

3.3. Mass Metrics and the Kinetic Energy. The potential energy U{Q) of 
([7]) is particularly simple, but the kinetic energy K{P) seems less natural. In this 
section we will see that it is related by duality to a Hermitian metric which will 
play an important role later on. 

Define a Hermitian mass metric on C'^ by 

(14) (y, W)^ — (mimat^i^ 14^12 + m^imiV^^W^n + ^27713^23^^23) • 

777 

The corresponding norm is given by 
m 

The mass norm 



(15) 



777l7772|Ql2p + 77737771 1 Q31 ^ + 7)727773 | (923 I 



r-= 101 = VWQ) 

provides a natural measure of the size of a configuration Q — {Q12, Q31, Q23) G C'^. 
In particular, r = represent triple collision. There is a dual m,ass metric on C"^* 
given by 



(16) 

with squared norm 
(17) 



(P, P) = 777 



Pl2^i2 
777 1 7772 



77737771 



^23-^23 
777270,3 



\P' = 



121 



\P:- 



311 



IP 



23 1 



777 1 777 2 7773 7771 



Note: Altogether we have three interpretations of 



77727773^ 

(., .) depending on whether the 



arguments are two vectors (eq 14), two covectors (eq 16), or a vector and a covec- 



tor, (eq 10). All three pairings are Hermitian, begin complex-linear in the second 



argument and anti-linear in the first. 

Introduce the notation Wo = yV\0 (and a similar notation for any vectorspace). 
If Q G Wo then it is easy to check that the vectors Q, N, T form a Hermitian- 
orthogonal complex basis for TqC"^ with respect to the Hermitian mass metric. 
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where 

Q = {Qi2, Q31, Q23) 

^-^g^ N = {1713,1712, mi) 

rp _ / Q3I _ _ Q12 Q12 _ Q3I N 

TO2 mi ' mi ma ' ma m2 

Q is a radial vector and N, T are respectively normal and tangent to W. Clearly 
{Q,T} is a basis for W. 

The next lemma shows the relationship between the kinetic energy and the dual 
of the mass metric. 

Remark. On terminology. A nondegeneratc quadratic form on a vector space, or 
on the fibers of a vector bundle, determines uniquely a quadratic form on the dual 
vector space, or on the fibers of the dual vector bundle. We refer to this dual 
quadratic form as either the 'cometric' or the 'dual norm'. 

Lemma 1. The kinetic energy satisfies 

where \P\ is the dual mass norm and where ttw : C'^ — > C"^ is orthogonal projection 
onto W with respect to the mass metric. 

Moreover, K{P) can be characterized as one-half of the unique translation- 
invariant quadratic form on TqC^ which represent the dual of the restriction of 
the mass norm to TqW. 

Proof. A direct computation shows that 

2 _ mm? ^ IA2-P31P IP23-P12P \P31-P23? ^ 

' ' |iV|2 2mi 2m2 2m3 ^ 

On the other hand, dual norms, or 'cometrics', can be characterized by the property 
that for any orthogonal basis {Q,N,T}, 

ipi2 ^ mQ)? , mm? , \{p,t)? 

' ' |Q|2 ^ |iv|2 |r|2 • 

Hence 



oK(P^-\P\2 m^m? _ \{p,Q)? , \{p,T)? 

ZAi^i-j — — -I- 



|iV|2 \Q\^ 

and this is also the formula for \P o ttw?- 

If we view TqW as the quotient space of TqC^ under momentum translations, 
then any norm on TqW is represented by a unique translation-invariant quadratic 
form on TqC^. In particular, this applies to the dual norm of the restriction of the 
mass norm to TqW. Since {Q, T} is an orthogonal basis for TqW with respect to 
the mass metric, this "lift" of the dual norm will be given by 

m^Q)\\\SE^ = 2KiP). 



QED 
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3.4. Parametrizing W. Let 

ei = (ai2, 031, 023) 62 = (foi2, 631, ^23) G W 
be any complex basis for W. The corresponding coordinate map is / : — > W C 

/(Ci,6) = 6 ei +6e2 or Q^j = Ci + 6 &y 
whiere ^ = (Ci;^2) G are the new coordinates. 

Extend / to a map F : T*C'^ ^ W x C^* by letting P e C^* be any solution to 
the equations 

(P, ei) = fji {P, 62) = ^2 
where = (771, 772) G C^* is the dual momentum to ^ and N is the normal vector to 
W from ( |18[ ). Any two solutions will differ by a momentum translation which will 
not affect the computations below. This definition makes F partially symplectic, 
where the symplectic structure on T*C^ derives from the canonical one- form 

= re(77,0 = re(?7i ^ +7726)- 

To find the new Hamiltonian, note that the pull-back of the Hermitian mass 
metric is 

Clearly this can be viewed as the pull-back of the restriction of the mass metric to 
W. The dual of this metric is 



9 



g = det G. 



.922 ^.921 
-512 511 . 

It follows from lemma[l]and the fact that the momenta also transform as pull-backs 
that the kinetic energy will be one-half of the dual norm. 
The Hamiltonian becomes 

(20) m.v)^\fG-'v-U{0 
where 

TO17712 TO1TO3 , '7l2m3 C ^ h C \ 

P12 P31 P23 

Example 1 (Heliocentric Coordinates). One can form such a parametrization of 
W by choosing one of the masses, say mi, to play the role of the origin. Set 

Q12 = —£,1 Q31 = £2 Q23 = £.1 — (2 

so that ^1,^2 G C are the coordinates of 7712,7713 relative to 7711. The corresponding 
basis for W 

ei = (-1,0,1) 62 = (0,1,-1) 
and the momenta fji = (P, e;) satisfy 

Vi — ^23 ^ P12 V2 = ^31 — -P23- 
For example, we can choose 

P12 = -m P31 = m P23 = 0. 
Substituting into Hred gives the familiar Hamiltonian 

Tjlc \ W+m? , \m? , \V2? 777 1 7772 777 i 7773 777 2 7773 



27771 27772 27^3 161 161 16 - 6 



PLANAR THREE-BODY PROBLEM 



13 



Example 2 (Jacobi Coordinates). Alternatively one can introduce Jacobi coordi- 
nates ^1, ^2 by setting 

Tn ■ 

Ql2 = -Cl Qsi = £.2 + i^26 Q23 = -6 + i^lCl 1^1 = 



nil + ra2 



This corresponds to the orthogonal basis 

ei = (-1, V2, Vl) 62 = (0, 1,-1) 

and we have mass metric 
Ml 



G = 



TOim2 {rni+m2)m^ 

wlicre /ii = /i2 = • 

mi + 7712 rn 





The momenta satisfy 

?7i = -Pv2 + V2P31 + i'iP23 m = ^^31 - P23 
and for an inverse we could choose 

P12 =0 /31 = ??i + i^im P23 = m- V2V2 



From ( 20 1 we get the equally familiar Hamiltonian 



Hi^ rj) - ^^^^^ + _ mim2 _ mimz ^2^3 

2^1 2/i2 ICll 16 + 1^261 \S.2-'^i£.i\' 

4. Spherical-Homogeneous coordinates 

The Hamiltonian Hrei{Q,P) of ([t]), representing the translation-reduced planar 
three-body problem, has further symmetries. The potential function U (Q) is sym- 
metric under simultaneous rotation of the Qij in C and is also homogeneous of 
degree —1 with respect to scaling. In this section we exploit the scaling symmetry 
by converting the system to spherical coordinates. This will be useful later when 
we blow-up the triple collision singularity. 

We use the mass norm r — \Q\ as a measure of the size of a configuration 
Q = {Qi2,Q3i,Q23) & C^. In particular, r = represent triple collision. For 
Q G Cq we want to define a spherical variable X G to describe the normalized 
configuration. However, instead of using the unit sphere = {X G : \X\ = 1} 
we will view the sphere as the quotient space of Cg under scaling by positive real 
numbers. This gives a convenient way to work globally on S^. We will take a 
similar approach when working with the complex projective space CP^ in the next 
section. 

Let M = T*Cg ~ eg X C^* with the standard symplectic structure. Let G = M+ 
be the group of positive real numbers and let G act on M by A: • {X, Y) — {kX, Yj k) 
where X e Cg, F £ C^*, fc > 0. We will use the notation [X],[X,Y] to denote 
equivalence classes under scaling. In other words, two vectors X, X' G Cg are 
equivalent, denoted X' - X if X' = kX for some k > 0. Similarly {X', Y') - (X, Y) 
if X' = kX, Y' = Y/k for some fc> 0. 

The momentum map for this group action is given by S{X, Y) — re(y, X) where 



the angle bracket denotes the Hermitian evaluation pairing ( 10 1. Fixing this scaling- 
momentum to be re(y, X) = and passing to the quotient space we get a reduced 
symplectic manifold which can be identified with the cotangent bundle T*S^. This 
is a special case of cotangent bundle reduction at zero momentum, as described in 
theorem [21 Introduce the notation 

T* C3 = S-\0) = {{X,Y) e T*Cl : re(y,X) = 0}. 
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Then we have 

We are going to pass from the relative configuration variable Q € Cq to a size 
variable r and a homogeneous variable X G Cq. 

Definition 2. {r,X) e x Cq are spherical-homogeneous coordinates for the 
configuration Q G Cg provided r = \Q\ and [X] = [Q]. 

X will be defined only up to a positive real factor and will be viewed as repre- 
senting a point of S^. We can use Q itself as a homogeneous representative of the 
corresponding point in S^. Hence we define a spherical-homogeneous coordinate 
map 

f:Cl-^R+xCl r=\Q\,X = Q. 
Extend f{Q) to a map F{Q, P), F : T*Cl T*R+ x T^^^C^ by setting 

\Q\ IQP ^ 

Here Pr gR.*,Y G €?* are the conjugate momentum variables to r, X and Q* is 
the dual covector to Q with respect to the mass metric. By definition, this means 
the unique covector in C^* such that 

{Q\V) = {Q,V) 

where the first angle bracket is the evaluation pairing and the second is the mass 
metric. We find 

(21) Q* = —(mim2Qi2,mim3(53i, 7712/713(323) e C^*. 

m 

A pseudo-inverse G{r,pr,X,Y), G : T*M+ x T*^^^^? T*C§ to F is given by 

(22) G: Q.r^ + 

We have G o F = id and 

F oG{r,pr,X,Y) = {r,pr,kX,Y/k) where fc = 7^. 

Hence foG = idmodM+. 

To check that F, G are partially symplectic, compute the pull-backs of the canon- 
ical one-forms 



(23) 9=prdr + re(yi2 dXi2 + Y31 dX^i + Y23 dX; 



23; 



and from (12 1. We find G*9 = Q while F*Q = 6 + . . . where the omitted terms 
are divisible by ve{Y,X). Hence the maps preserve the restricted symplectic forms 
as required. 

The spherical-homogeneous Hamiltonian is Hgph — Hrei ° G. Using the formula 
for Q in (22), the potential U{Q) becomes Usph{r,X) = ^V{X) where 

/mi 7712 77737/7,1 777277t3' 
1, IX12I ^ IX31I ^^2^ 



(24) V{X) = \X\U{X) = \X\ 



Note that V is invariant with respect to scaling of X so it determines a well-defined 
function, F : — >■ M, which we will sometimes write as y([X]). 

The kinetic energy is Kspk — K{P) where P is given by (22 1. It follows from 
lemma [l] that the two terms in ( 22 ) are orthogonal with respect to the quadratic 
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form K. To see this, note that they are orthogonal with respect to the dual mass 
metric since {Y,X*) = {Y,X) = 0. Since X e W we have 

(y O TTW, X* O TTw) ^{Yo TTW, T^W X) = (F, X) = Q 

SO X* o TTw, 5^ o Ti'w are still orthogonal. Evaluating K separately on the two terms 



of (22) we find 



(25) K,,h - \pI + ^-^K{Y) 
and so the spherical-homogeneous Hamiltonian is 

(26) Hsph{T,Pr,X,Y) = + ^if(r) - IviiX]). 

Theorem 5. The Hamiltonian flow of Hsph on r*M+ x T*Cq has invariant sub- 
manifold {re(y,X) = 0} and the quotient of the restricted flow by the scaling sym- 
metry is equivalent to the Hamiltonian flow of Hj-ei on r*CQ. This submanifold con- 
tains a codimension 2 invariant submanifold {re(y, X) — 0, X12 + X^i -\- X23 — 0} 
for which the quotient of the restricted flow by the symmetry of scaling and trans- 
lations of the Yij is conjugate to the flow of the zero total momentum three-body 
problem reduced by translations. 

Proof. For the first part we apply theorem [l] with Ah = T*C§, M2 = r*M+ x T*Cg 
and symmetry groups Gi = {id} and G2 — K^- The momentum level is S{X, Y) = 
re(Y,X) = 0. It was shown above that the maps F, G between T*Co and 5*^^(0) 
are partially symplectic pseudo-inverses. 

For the second part we change the groups to be Gi = C* and G2 is a semidirect 
product of the scaling group M"*" and the momentum translation group C* with 
group multiplication 

(^2, C2) • (fci, Ci) = (fc2fcl, Ci/fc2 + C2) 

where {ki,Ci) £ 1R+ x C*. The momentum levels are {Qtot — 0} and {Xtot ~ 
0, re(y, X) = 0} respectively, and these are fixed by the actions of the groups. The 
maps F, G restrict to maps between these level sets and the restrictions are partially 
symplectic pseudo-inverses. QED 



If we use the formula K{Y) = ^Y^ BY with B from (|9j) wc find Hamilton's 
equations for Hgp^ are 

r ^Pr 

2\X\^K{Y) 1 



,V{X) 



(27) . |XP 

X = ^—}-BY 



Y ^ \DViX) - 

The quotient space of r*K+ x T*p^Co mentioned in theoremjsjis diffeomorphic to 
T*R+ X T*S^ (by simply thinking of X, Y as homogeneous coordinates for [X, Y] e 
r*S5). The quotient space of r*M+ xT/^^^Ci^ is diffeomorphic to T*m+ y.T* S{yV) 
where SiW) = W n is diffeomorphic to S'^. Hence the reduced space is eight- 
dimensional as before. The reduced flow is just the translation-reduced three-body 
problem in spherical coordinates. 
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At this point, instead of reducing the number of dimensions, we have actually 
increased it from 12 to 14. The value of the present formulation lies in the fact 
that it has been put in a form where double collisions can be easily regularized 
and the triple collision easily blown-up without destroying the symmetry among 
the masses. As in the previous section, one could explicitly realize the reduction to 
eight dimensions by parametrizing the subspace W. However we will not do this 
here. 

5. Reduction by rotations - the shape sphere 

Next we form the quotient by rotations. Since we are using complex coordinates, 
the combined action of scaling Q by a real factor r > and rotating Q by an angle 
9 is represented by Q h- > kQ where k = re^^ G Co = C \ 0, the space of nonzero 
complex number. A point in the resulting quotient space represents the size and 
shape of a configuration. 

5.1. Radial-homogeneous coordinates. As before we will measure the size by 

r = \Q\. 

To represent the shape we project Q G Cg to the quotient of Cg by the action of Cq. 
This quotient space is the complex projective plane P(C'^) = CP^. Homogeneous 
coordinates will provide a way to work globally on the projective plane, just as 
they did for the sphere in the last section. For X & Cl let [X] G CP^ denote 
the corresponding element of the projective plane, i.e., the equivalence class of X 
under the relation that X ~ Q if X = kQ for some A; G C, 7^ 0. (Thus the square 
bracket will now mean a projective point rather than a spherical one.) 

Definition 3. (r, X) are a pair of radial-homogeneous coordinates for Q G Cq if 
r= 101 and [X] = [Q] G CP^ 

X is defined only up to a nonzero complex factor. We can take X = Q itself to 
define the radial-homogeneous coordinate map 

f:Cl^R+xC^ r=\Q\,X = Q. 

Remark. Despite the fact that spherical-homogeneous coordinates and radial ho- 
mogeneous coordinates are both denoted {r,X) there are differences between the 
two coordinate systems. Spherical- homogeneous coordinates represent points in 
Cq ~ R+ X whereas radial homogeneous coordinates represent points in the 
quotient space {CD/S'^ ~ M+ x CP^ 

If we include the origin and form the quotient space under rotations we have 
C^/S^ = Cone(CP^), the cone over CP^, where the cone point corresponds to total 
collision G C'^. For any topological space X, we can form the space Cone{X) 
which has a distinguished cone point * and Cone{X) \ * = M+ x X. In this case, 
the cone is not a smooth manifold. 

The equivalence class [X] — [Q] G CP^ represents the shape of a three-body 
configuration only if Q G W. Restricting to such Q we get [Q] G P(W), where 
P(>V) is the projective space of the subspace W C C^. Since W is a two-dimensional 
complex subspace, P(H') is a projective line, i.e., 

P(W) ~ CP^ ~ S^. 

P(W) will be called the shape sphere. 
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Any function on our original configuration space which is invariant under trans- 
lation, rotation, and scaling induces a function on the shape sphere, the most 
important example being our homogenized potential 

V{X) = \X\U{X) : PW ^ M. 

We will also use homogeneous momentum variables. A pair {X,Y) G T*Cq ~ 
Cq X C^* win represent a point of T*CP^. Let G = Co be the group of nonzero 
complex numbers and let G act on T*Cq by k ■ {X,Y) — {kX,Y/k). We will use 
the notation [X, Y] to denote equivalence classes under scaling. In other words, 
{X', Y') - {X, Y) if X' = kX, Y' = Y/k for some nonzero keC. The momentum 
map for this group action is given by the Hermitian evaluation pairing a{X, Y) = 
{Y,X) E C. The real part of the complex number a{X,Y) is the real scaling- 
momentum S{X, Y) (which we want to be zero as in the last section). On the other 



hand, from ( 11 ) we see that ima{X, Y) — —i /i, where /i is the angular momentum. 

If we fix the complex scaling-momentum to be (F, X) — and pass to the 
quotient space, then as in theorem [2] we get a reduced symplectic manifold which is 
naturally identified with the cotangent bundle T*CP^ with its natural symplectic 
structure. Introduce the notation 

T;X' = a- 1(0) = {(X,Y) eT*Cl: {Y,X) = 0}. 



Then we have 



t;^C7Co ~ T* 



If, on the other hand, we fix the complex scaling-momentum to be (Y, X) = —i fj, 
and pass to the quotient space we still get a reduced symplectic manifold which 
can be identified with the cotangent bundle T*CP^ but with a twisted symplectic 
structure, as described in theorem [3] More about this below. 

To get a system equivalent to the reduced three-body problem we will also need to 
include the radial variables. Restrict X to W and quotient by the action of the group 
C of translations in F-momentum space. Let M = T*M+ x T*Cq with coordinates 
{r,Pr,X,Y) and let G = Cq x C acting by {k,c) ■ {r,pr, X,Y) ^ {r,pr,kX,c- (Y/k)) 
where c ■ Y — (Y12 + c, Y31 + c. Y23 + c). Fixing the momentum level J{X, Y) — 
{u{X, Y),Xtot) = {—i M, 0) G C'^ and passing to the quotient space gives the reduced 
phase space 

P^{ir,pr,X,Y) : = -^M,Xl2 +^31+^23 = 0}/G 

of real dimension dim P = 14 — 4 — 4 = 6 as expected. In fact we have 

P ~ T*M+ X r*P(W) ~ T*M+ X T*S^. 

We still need to find the reduced Hamiltonian and show that the reduced Hamil- 
tonian system is equivalent to the reduced three-body problem. This is easy to do 
starting from the spherical Hamiltonian in the last section. Indeed, the passage 
from the homogeneous-spherical variables {r,pr, X,Y) G T*M+ x r*CQ to the cor- 
responding radial-homogeneous ones is just given by the identity map. The new 
feature here is that the symmetry group is enlarged from K"*" x C* ~ K"*" x C to 
Co X C. Then we have the following extension of theorem [Sj 
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Theorem 6. The Hamiltonian flow of Hgph on r*M+ x T*Cg has an invariant 
set where (Y.X) = —ifi. The quotient of the restricted flow by the complex scaling 
symmetry is equivalent to the Hamiltonian flow of Hrei on r*CQ/S^. There is 
another invariant set where (Y, X) — —i fi and ^12+^31+^23 = and the quotient 
of the restricted flow by the complex scaling symmetry and by translations of the Yij 
is conjugate to the flow of the three-body problem with zero total momentum and 
angular momentum fj,, reduced by translations and rotations. 

Proof. The maps F and G as in the proof of theorem [5] restrict to maps of the 
/i angular momentum levels. They are still partially symplectic pseudo-inverses. 

QED 

The next step is to use a momentum shift map to pull-back the problem to the 
zero-angular-momentum level. This expresses all of the reduced problems on the 
same phase space and makes the role of the angular momentum constant explicit. 
Let 

(28) <!>^{r,pr,X,Z) = ir,pr,X,Y) Y = Z + ^,T{X) T{X) = '^* 



where 



1 



im{Y,X) = im{ifi—^,X) = -^re(— = -fi. 



X* = —(mim2Xi2, msmiXsi, 77127713X23) G C''*. 
m 

Note that : J-^{0,0) J-i(-i^,0) since if {Z,X) = we have 

Composing H^ph with $p we get a Hamiltonian 

(29) H,ir,pr,X,Z) = ^ipl + ^) + ^-^KiZ) iy([X]). 

To verify this we need to show that the kinetic energy can be written 



(30) K,^-ipl+'^)+^-^K{Z) 



This decomposition follows from an orthogonality argument based on lemma [T] 
Namely, the vectors i X and Z are orthogonal with respect to the mass metric and 
the first one lies in W. Then, as in the last section, lemma [T] shows that they are 
orthogonal with respect to the quadratic form K and so K{Y) = K{pr{X))+K{Z). 
K{iJ,T{X)) gives ^^-term in Kf^. 

Equation [30] gives a decomposition of the kinetic energy into radial and angular 
parts and a third term which can be viewed as the kinetic energy due to changes 
in the shape of the configuration. Some authors call this decomposition of kinetic 
energy, or the consequent orthogonal decomposition of velocities the "Saari decom- 
position". (See j25].) In the next subsection we show how this last shape term can 
be understood in terms of the Fubini-Study metric on the shape sphere. 

5.2. Fubini-Study Metrics and the Shape Kinetic Energy. Using a complex 
orthogonal basis, we give a simple decomposition of the dual mass metric which 



leads to deeper insights into the kinetic energy decomposition (30). Since the 
shape sphere has complex dimension one, there are some very simple formulas for 
the shape term of this decomposition. 
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To describe the Fubini-Study metric (also called the Kahler metric), let V denote 
any complex vector space and let (V, W) be any Hermitian metric on V. If X G 
Vo = V \ then the corresponding Fubini-Study metric on TxV is: 

(31) {V,W)^s = . 

As a bilinear form on TxV, the Fubini-Study "metric" is degenerate with kernel 
the complex line spanned by the vector X. But it induces a bona fide Hermitian 
metric on the projective space P(V). 

To see this, let tt : Vq — >■ P(V) denote the projection map: ii{X) ~ [X]. The 
tangent map Ttt : TVq rP(V), Ttt{X,V) = {[X], Dtt{X)V) has the property 
that Tii{X, V) = Tti{X\ V) if and only if X' = kX and V' ^ kV + IX for some 
complex numbers k 0,1. So it is natural to view the tangent bundle TP(V) as the 
set of equivalence classes of [X, V] of pairs {X, V) E Vq x V under this equivalence 
relation. It is easy to check that the formula for {.,.)ps is invariant under this 
equivalence relation and so it gives a well-defined Hermitian metric on P(V). The 
real part re(y, W) ps gives a Riemannian metric on P(yV) and the imaginary part 
gives a two-form called the Fubini-Study form which will be important later 

nFsiV,W)^im{V,W)FS- 

Starting with the mass metric on V = C"^ we get a Fubini-Study metric on 
CP^. However, because of lemma [l] we will be interested in its restriction to the 
two-dimensional complex subspace W C which we denote by (., ■)fs,w, which 
induces a Hermitian metric on the shape sphere P(W). 

Our goal is to show that the shape kinetic energy is the cometric dual to this 
Fubini-Study metric on P(yV). (By a 'cometric' on a manifold X we mean the 
fiberwise quadratic form on T*X which is dual to a Riemannian metric on X.) To 
this end we will need to describe cometrics on projective space in homogeneous 
coordinates. We continue to identify r*CP^ with the quotient space of T*j.C'^ = 
{{X, Z) G CqX C^* : {Z, X) — 0} under the complex scaling symmetry. In the same 
spirit, the cotangent bundle T*V{yV) is the quotient space ( a symplectic reduced 
space) 

r*p(w) ~ {t;^^^cI)/Co x c 

where 

T;r,wCl = {(X, Z) e W X C^* : (Z, X) = 0, X ^ 0} 

and where the group Cq x C represents the scaling symmetry and the momentum 
translation in Z-space. We refer to {X,Z) as homogeneous coordinates on P(>V). 
The restriction of Z e C^* to W representing a covector in Tj*^jP(yV). Expressed 
in homogeneous coordinates a cometric on P(W) is a function of the form Q{X, Z) 
which is quadratic in Z and invariant under the Cq x C action. 

Theorem 7. The Fubini-Study cometric \ZWgYV i-^] ^ related to the 

kinetic energy (formula {Kf^ ) by 

1 

Proof. Substitute (X, Z) for (Q, P) in formula Use {Z, X) = to get K{Z) = 
^l^|2 {Z, T) . The vector field T{X) appearing in that formula is tangent to W 



l^lkw - \X?K{Z) 
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and orthogonal to X, hence fits the hypothesis of lemma [2] immediately below. The 
lemma asserts that with e{X) = |J^^T(X) we have |^||s.w = K^' '^{X))\^ ■ QED 

Lemma 2. Let T{X),X G Wo be a nonzero complex vectorfield tangent to Wo 
and normal to X with respect to the Hermitian metric mass metric. Then e{X) = 
I X I 

|7jip^T(X) is a unit tangent vectorfield on Wo with respect to the pulled back 
Fuhini-Study metric {■,.)fs.W- Moreover 

(32) {V, W)fs.w = ^' ^ ^ W/(CX) - Tyxfm. 

and the pulled-back cometric is given by the quadratic form 

(33) \Z\%s,w = \{Z,e{X))\' Z e T^^^.C^. 



Proof. Since T{X) is orthogonal to X, (31) gives 

and so e{X) is a Fubini-Study unit vector at X. 

The tangent space T^W has complex dimension two and {X,e{X)} is a basis. 
If we expand V e TxW a.s 

^- W WIT™ 

and similarly for W, then since X is in the kernel of (., .) ps we get 

{V, T{X)){T{X),W) {V, e{X)){e{X),W) 



{V,W)fs,w^{V,W) 



FS 



|X|2|r(X)|2 |X|4 

as claimed. 

Observe that if E, (•, •) is a one-dimensional complex Hermitian vector space with 
unit vector e then the cometric on E* is given by the quadratic form Z e E* i-> 
\{Z, e)|2. From this observation the last formula of the lemma follows. QED 

Remark. The manifold P(W), being a two-sphere, admits no non- vanishing vector 
field. So how did we just construct a unit vector field e{X) to this two-sphere? 
We did not! The gadget e{X) is a unit section of the pull-back /*rP(W) of this 
tangent bundle by the homogenization map / : Wo — >■ P( W) which sends X [X] . 
This pull-back bundle can be viewed as a sub-bundle of TWo, and hence e{X) is a 
vector field on Wq. 



Using the vector field T{X) of formula ( 19 1 (with X substituted for Q, we obtain 
the Fubini-Study unit tangent vector 



mim2TO3 


fX^i 


X23 


X23 


X12 


X12 


X31 


m 


V m2 


nil 


mi 


ma 


ma 


J7l2 



e{X) = 

From this expression we get simple formulas for the Fubini-Study metric and two- 
form on W: 

.X , > mimsma _ mim2ma . _ 
(34) {.,.)fs,w^ |y|4 cr(8)cr Hfs.w = ima'^a 

where the complex- valued one-form a is given by any of the following formulas 
(35) 

<T = (e, dX) = Xai(iXi2 — Xi2dX3i — Xi2(iX2a — X2a'iXi2 = X2a'iX3i — Xai(iX23. 
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For example, the second formula for a is obtained by eliminating X23,dX23 from 
{e,dX) using the equations X23, = —^12 ^ -^31 and dX2z — —dXi2 — dX2,i. Note 
that the formulas for a are independent of the masses. This implies that the Fubini- 
Study metrics for different masses are all conformal to one another. 

Similarly we get a formula for the dual norm and the shape kinetic energy: 

(36) \X?KiZ)^\\Z\ls^^ = p^ 

where a{Z^ is given by any of the following formulas 

(37) ^ 

a = — (TOirn2-'^12(^23 - ■Z'3l) + "l3"llX3i(Zi2 - ^23) + ™2TO3-'^23(^31 - ^12)) 

m 

|X|2(Z31 - Z12) |X|2(Zi2-Z23) |X|2(Z23 - Z31) 



-''^23 X12 

Our identification of the shape kinetic energy with the Fubini-Study cometric 
gives an alternative formula for the reduced Hamiltonian on T*^C^ 



(38) H^ir,Pr,X, Z) = \{pI + ^) + ^\Z\ls^y, - ^-V{X). 

where \Z\'^pgy^ is the Fubini-Study cometric on W. 



5.3. Induced symplectic structure and the reduced differential equations. 

Using the momentum shift map, we have pulled back the Hamiltonian to the re- 
duced Hamiltonian defined on the zero-angular momentum level r*M+ x T*^£,^ 
where 

TpVC^ = {{X,Z) £ T*C^ : {Z,X) = 0}. 

However, as described in theorem |3] there is also an induced symplectic structure 
on this set which different from the restriction of the standard one. The pull-back 
of the canonical one form 6 under the momentum shift map ( 28 ) is 

$* 6* ^prdr + ie{Z, dX) + im{X* ,dX)=(d + ^Qi 



with 



. {X*,dX) . {X,dX) 
Bi = im — — -pr — = im ■ 



where we changed the evaluation pairing to the mass metric in the second equation. 
The modified symplectic form will be 

fi^ — — fidQi 

where we find 

.^q^ ,n 9- {dX,dX)\X\^~{dX,X){X,dX) 

(39) dei = 2im 1^^^ ^2nps 

where i^'pg is the Fubini-Study two-form determined by the mass metric on C'^ (as 
opposed to its restriction to W as in section 5.2 Geometrically, ^'pg represent the 
curvature of the circle bundle — CP^. 

Once we have fi^ we calculate Hamilton's differential equations using the defining 
equation for Hamiltonian vectorfields: 

(40) ir,pr,X,Z)^n^^dH^. 
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The interior product with the standard form gives the usual result: 

{r,Pr,X,Z)_in = -pr dr + r dpr - ie{Z,dX) + re{X,dZ). 

Since i^'pg involves only dX, it can be viewed as a two-form on C'^ instead of on 
phase space. Moreover, it only affects the differential equations for Z. Hamilton's 
equations read: 



r 


— ^M,Pr 


Pr 




X 


= Hf^^z 


Z 





(41) 

^ — ^^n,z 



where i?^ is given by ( 29 ) . The term involving the Fubini-Study metric will be 
called the curvature term, T'^^^^ = —2^,H^ z-^^'fs- 

Lemma 3. // X e W and {Z, X) — 0, then the vector H^^z is in W and 
{X,H^^z)^O.Infact 



(42) H,,^z - e e W 

where e[X) is as in lemma^ 

The curvature term T^^^^ is equivalent under the translation symmetry in C^* 

to 

(43) Tcurv = -^iZ. 



Proof. From (29 1 we have = ^-^DK{Z). Note that since Z e C^* , we have 

DK{Z) : C^* M. By duality we can view DK{Z) as a vector in C^. Let X e W. 
Since X — H^ z and W is invariant, we must have H^ z <= W. If {Z,X) = 
then an orthogonality argument as above shows K{Z + X*) = K{Z) + K{X*) 
which implies, since if is a quadratic form, that DK{Z){X*) ~ {DK{Z), X) = 
as required. 

In section 5.2 we showed that in the subspace {Z : {Z,X) = 0} we have 
\X\'^K{Z) = ^\{Z, e)p. In fact, we will see that the Z-derivatives of these two 
functions also agree: 

(44) \X\'^DK{Z) ^ (Z,e)e. 



To see that (44) indeed holds, note that differentiation along the subspace shows 
that they must agree when evaluated on any SZ with {SZ, X) — 0. On the other 
hand, both sides vanish on the complementary covector Z' = X* . Note that the 
right hand side was calculated, as always, by converting to real variables, finding 
the real derivative and then converting back to a complex vector. Equivalently, we 
expand 

^\{Z + dZ,e)f^^\{Z,e)f+ re{SZ, Jz;^) e) + . . . 

for all 5Z , showing that the vector in question is the complex representative of the 
real vector derivative. 

To show the equivalence of T'^^^^ and T^urv we will show that they agree when 
restricted to W. The argument can be based on a kind of Fubini-Study duality. 
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Namely, if G W we will show that 

(45) {H^^z,V)fs^^{Z,V) 

which means that r^H^ z is a dual vector to Z with respect to the Fubini-Stiidy 
metric on W. To see this note that (44) gives 



1 ((Z,e)e,T/) _ {Z,e){e,V) 
{H,.z,V)fs - - ,2|;^|2 ■ 

On the other hand any e W is linear combination 

Since e is a Fubini-Study unit vector, we have |e| — \X\ and so 

' ' r2|eP r2|X|2 



and (45) holds. From this we can calculate that for any y g W 

r^^Ci^) = -2Ai im(i/^,z, y)j^S = im(Z, y) - re(i Z, 1/). 

This means that T'^^rv T^urv agree as real-valued one-forms on W as claimed. 
Replacing T'^^^^ by Tcurv introduces only an irrelevant translation of the momentum 
Z. QED 

Taking this lemma into account we finally get Hamilton's equations for the re- 
duced Hamiltonian in the form 

r ^Pr 

. ^^' + \X\■'2K{Z) 1 

Pr = 1 2^\^) 

Applying theorem [l] to the momentum shift map and remembering theorem [6] 
we have: 

Theorem 8. The Hamiltonian flow of i?^ on T*R+ X T*Cl has an invariant set 
r*M+ X T*^€? where {Z, X) — with symplectic structure given by the restriction of 
the standard form minus 2iJ,Qps. The quotient of the restricted flow by the complex 



scaling symmetry is equivalent to the Hamiltonian flow of H on T*Cg/S^. Th 



ere is 



another invariant set r*M+ x T*^ yyC'^ where {Z, X) = Q and X12 + X^i + X23 — 
and the quotient of the restricted flow by the complex scaling symmetry and by 
translations of the Zij is conjugate to the flow of the three-body problem with zero 
total momentum and angular momentum p,, reduced by translations and rotations. 

This Hamiltonian system represents the reduced three-body problem in a way 
which is convenient for regularization of binary collisions and blow-up of triple 
collision. However, the phase space is still 14- dimensional. Next we describe how 
to find lower-dimensional representations of the reduced three-body problem by 
parametrizing the shape sphere in various ways. 
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5.4. Parametrizing the Shape Sphere. The shape sphere is the projective space 
P(yV). As in section 3.4 choosing a complex basis {61,62} for W gives a map 
/ : — W, X = /(C). By viewing X € W and ^ € as homogeneous coordinates 
we get an induced parametrization of the shape sphere fpr : CP^ — )■ P(>V). 

The formulas of section 3.4 (with (Q, P) replaced by {X, Z)) allow us to find the 
reduced Hamiltonian for any such basis. If 

61 = (ai2, asi, 023) 62 = (foi2, 631, ^23) e W 
then we have, as before, 

Xij — £,1 aij + £2 bij 

and 

77i = (y,6i) ^2 = (y",62). 

We define a Hermitian mass metric and dual mass metric for C, 77 to be the pull-backs 
of the metrics for X, Y. The squared norms are 



V 



where G is the matrix with entries Gij = {ci, ej) , and these squared norms represent 
the mass metric and cometric on W. 

The relation between the cometric and kinetic energy yields the Hamiltonian. 
(See equations ( 29|30| ) and theorem[7|) 

(47) H,ir,p^,£,,) = l{pl + ^ + ^) - IviO 

where the shape potential is 

V P12 P31 P23 / 

To make the map F of section [3^ be partially symplectic we need to alter the 
standard symplectic form in (C, 77)-space by subtracting 2/1 F*Q'pg. Pulling back 
the Fubini-Study metric (., .)fs by / gives the Fubini-Study metric in ^ space 



With the help of ( 34 1 one can show 

(■, ■)fs = 1^ CTo <^ ctq where ctq = £id£,2 - g = det G 

The Fubini-Study two-form is the imaginary part. 

Since cto is independent of the choice of basis, the Fubini-Study metrics for vari- 
ous choices of basis are all conformal to one another. If we choose an orthonormal 
basis the metrics are Euclidean. The Fubini-Study metric for a general basis is 
related to the Euclidean one by 



where the conformal factor is 
(48) KiO = 

where |ClLc = ICiP + 



FS = l^iO ^FS.euc 

9\£.\t 



I euc 
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The curvature term can be calculated directly from the definition H^^ri-^^FS 
and we find 

T _ 2/i . 

curv — "^V' 

Hamilton's equations in T*M+ x T*^C'^ are 

There are still 10 variables but the invariant set T*K+ x T^^C^ with {rj, 0=0 
is 8-dimensional and we have a complex scaling symmetry. The introduction of an 
affine coordinate on the projective line yields a full local reduction to 6 variables. 
For example, consider those points [^] = [^1,^2] G CP^ with ^1 ^ 0. If p is any 
nonzero constant complex number then every such point has a unique representative 
of the form 

[^1,^2] = [p,z] z = x + iyGC, 
thus parametrizing almost of the shape sphere by a single complex variable z, the 
affine coordinate. Of course the roles of ^1 , ^2 could be reversed to parametrize the 
subset with ^2 7^ 0. 

If C = a + 1/3 G C* denotes the momentum vector dual to z then the unique 
extension of f{z) = {p,z) to a partially symplectic map T*C — > r*^C^ = = 
0} is defined by 

6 = P ^2 = z r?i = -zC/p m = C- 
One computes the mass metric is 

1^(^)1^ = giAP? + 922\z\'^ +2Te{pgi2z) 

and the cometric is 

This gives a Hamiltonian system with 6 degrees of freedom 

(50) H^ir,pr,x,y,a,l3) = Upl + ^l+ '^^f '"'f ) - ^-V{x,y) 

where 

V{Z) = 1^(2)1 — + — + — p^J = \aij + bijZ\. 

\ P12 Psi P23 / 
The Fubini-Study form is 



^ 9 . ,_ , gdx A dy 
ilfrq = '- im dz ^ dz = . 

' \pmzw \pmzw 

The curvature term is just 

J- curv — ^S> 

as usual. 
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Example 3 (Projective Jacobi Coordinates). As a first example, consider using 
Jacobi coordinates as in section |3.4[ only this time applied to the homogeneous 
variables X, Z. As before, the basis which defines the Jacobi coordinates is the 
orthogonal basis 

ei = (-1, V2,v{) 62 = (0, 1,-1). 

We have 

^ = (-6:6 + V2i\,-i2 + Vxi\) i = (-Xi2, 1^1X31 - V2X2Z) 

and 

Z = {0,r]i + iyir]2,r]i - V2r]2) V = (-■Z'12 + 2^2.^31 + ;^i2'23, Z31 - Z23) 

where, as usual, Z is non-unique. 

The Hamiltonian is ( 47 ) where the shape potential is 

'r7iiTO2 miTOa ^27713 



^(0 = ICI 



161 16 + ^^261 16-^^161 

The mass matrix G = diag(/ii, /i2) has determinant g — fiiii2 = "^^"^"^^ and 
associated norm and conorm: 

ICP = ml6P + /.2|6P H^ = M + M. 

Ml M2 



Hamilton's equations with the curvature term are given by (49). 

If we introduce affine variables by setting S.i = p,^2 — z as above and if we 
choose p = .J— the mass norm reduces to 

|6' = M2(i + x2 + y') 

and we get the affine Jacobi Hamiltonian 

H^{r,pr,x,y,a,P) = - ( + — H — I - -V{x,y) 

Hamilton's equations with the curvature term are 
r ^Pr 



(1 + + ?;2)2 

X = ^ — '- a 

(51) (1 + x^+yYn 

y = P 

a = ^V^x, y) - 1(1 + x2 + y2)(^2 ^ ^2)^ ^ 
/j = IVyix, y) - 1(1 + + y2)(«2 ^ ^2)^ _ 

Example 4 (Equilateral Coordinates). In projective Jacobi coordinates (6j6): 
the binary collision points 612, 613, ^23 a-i'e located at the projective points 

[1,0], [1,-1^2], [i,i^i] ecpi 

while the equilateral triangle configurations (the Lagrange points) are at 

\i,e±] e cpi 
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where 



1711-1712 , V3 . 1^1 - 1^2 , V3 . 

t+ = — ^ ± — I ~ ± — I. 

2(toi+TO2) 2 2 2 



Using a Mobius transformation, we can put three points anywhere we like on the 
shape sphere, CP^. Remarkably, it turns out that if we put the binary collisions at 
the third roots of unity 

(52) [Ci,6] = [i,i],[i,^],[i,cs]ecpi 

1 

with Lu — —7^ + mf-, then the equilateral points are automatically moved to the 
north and south poles 

[1,0], [0,1]. 

These coordinates were introduced in [22 . 

These coordinates are obtained by choosing the basis 

ei^{l,uj,ui) 62 = -ei = (-l,-aj,-a;) 

for W. The coordinate change map is X = ei + ^2 ^2 or 

-^^12 = Cl ~ C2 -^31 = '^^l — X23 = wCl - ^^^2- 



and indeed takes the roots of unity (eq. 52 ) to the binary collisions. Setting 
^2 = we see that IX12I = IX32I = [-''^23 1 corresponding to an equilateral triangle, 
with the same result if ^1 =0. Thus the coordinate change map sends the poles 
^ = [1,0], [0, 1] to the equilateral triangles. 

The mutual distances (of the homogeneous variables) pij — \Xij\ which appear 
in the shape potential are very simple: 

Pl2 = |Cl-6l P31 = 1^1-^61 P23 - t^^l- 

The mass metric can also be written in terms of these 

[C[^ = ^(m-i"i2P?2 + msmipl^ + 77127713^23)- 
It is represented by the matrix G with entries gij — (ei, 62): 

77717772 + 77737711 + 77727773 _ 171x1712 + 77737771 LO + 17121713 UJ 

gil — 922 = 3l2 = 521 = 

777 777 

and determinant g = det G = Smi^iama . 
The inverse transformation is given by 

6 = ^(^12 + '^^31 + UJX23) 6 = ^^(^12 + W^31 + WX23) 

and the momenta satisfy 

rii — Z12 + i^Z^i + WZ23 772 = -Z12 - 0JZ31 — U1Z23- 



Choosing afhne variables by setting £,1 = z,^2 — ^ we get the Hamiltonian (50) 
with 

[C(2)[^ = (777 l7n2[z — 1[^ + 7n3777l[z — w[^ + 77727773[z — aj[^). 

m 

The complexity of mass norm is perhaps outweighed by the fact that the potential 
is given by the wonderful expression 



771,11712 777i7773 77727773 

|z — 1| [z — wl Iz — w| 
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The advantage of these coordinates is that they provide the homogenized poten- 
tial V with "radial monotonicity" ' . Let E = x^+y^he the radial vector field 
in the z plane, where z ^ x + iy. Then E[V] > for < \z\ < 1, E[V] < for 
\z\ > 0, and E[V] = if and only if |z| = 1 or z = 0. (See Proposition 4 of [12]) 
This monotonicity was the key ingredient to the main theorem of |23| . 

5.5. Making the Shape Sphere Round. Instead of using projective or local 
affine coordinates, one can map the shape sphere to the unit sphere in K^. First we 
do this homogeneously, then restrict to the unit sphere to get another version with 
6 degrees of freedom. Let ^ — {£,1,^2) G be coordinates associated with some 
choice of basis ei, 62 for W. 

Consider the Hopf map h : C"^ ^ M.^ given by 

wi = 2re^^2 W2 = 2im^^2 W3 = |fip-|6P 
Using the Euclidian metric for w we get 

It follows that 

We will need formulas for pij = \Xij\ = \aij£^i + 6ij^2| in the variables Wi. We 
have 
(53) 

pI = |a..fiaP + h.nU' + 2re(665.,&.,) 

= i(|ajj|^ + |foijp)|w| + li\aijf - \bijf)w3 + re{aijbij)wi - iT[i{aijbij)w2- 



Then the mass metric will be given by 

(54) = ^(miTO2Pi2 + msmipl^ + 77127713^23)- 

If we let ai, a2, 03 be dual momentum variables, we can extend the Hopf map h 
to a partially symplectic map F : T^^C^ — >■ T*pf^M.^ by defining its (pseudo) inverse: 

77 = a o Dh :~ Dh*a. 

To find the reduced Hamiltonian in w coordinates we will exploit the fact that 
Euclidean metric transforms nicely. Recall that the shape kinetic energy is the 
dual of the Fubini-Study metric and that the latter is related conformally to the 
Euclidian metric with conformal factor where k is given by (48 1. In other 
words, since we are restricting to (77,^) = we have 



One can verify that the Euclidean norms transform under the Hopf map in such a 
way that 

where we are using the Euclidean norm on , M^* . Hence the reduced Hamiltonian 
on the sphere is given by: 

H^ir,pr,w,a)^ -[pI + '-^+ ' ' ^ ) - -V{w) 
2 \ K\w)r'^ J r 
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where \w\'^ ^ wf+wl+wl and |ap = + + Q^i and where the shape potential 
is given by 

f mim2 mim2 mim2^ 



V{w) = \aw)\ 



V Pl2 Pl2 Pl2 



with the Pij and |^| as in ([53]) and (54 1 



The Fubini-Study form becomes a multiple k/4 of the Euclidean solid angle form 



K 



^FS — Ti — 7^{widw2 A dws + W2dwz A dwi + w^dwi A dw2)- 
4 w P 



This leads to the curvature term 



2p 

Tcurv = I — n^a X w 



where w x a denotes the cross product in . 
The differential equations are: 



r =Pr 



(55) 

w = —a 

■ 1^,., ^ 4|a|2 4|zi;|2|a|2 2/i 

a = -DV[w) H ^-^3 — K„ + 1 — -^a x w 

r Kr-^ K^r^ |w|r^ 

From theorem [l] if we restrict to r*E+ x T*^f^M? = {(a, w)e«c = 0} and quotient 
by the scaling action of K"*", we get a reduced system equivalent to the reduced three- 
body problem. But (a, w)euc — implies that |?z;| is constant under the flow. Hence 
we have a 6-dimensional invariant submanifold given by \w\ = 1, {a,w)euc — 
representing the reduced three-body problem. The reduced phase space is T*IR+ x 
T*S2 and the shape sphere is represented by the standard unit sphere. 

To get to 6-dimensions with no constraints one could parametrize the sphere with 
two variables. If this is done with stereographic projection, the result is the similar 
to the affine coordinate reduction of section 15.41 On the other hand one could also 
use spherical coordinates 9, (j). However, both of these are just local coordinates 
while the system above is global, albeit constrained. 

Example 5 (Jacobi coordinates on S^). If we choose an orthonormal basis for W 
then we get the conformal factor k — 1 and the resulting Hamiltonian will have a 
simpler shape kinetic energy. For example we could normalize the Jacobi basis of 
example [3] to 

^ :(-l,Z.2,^^l) e^, = ^(0,1,-1). 



The coordinates £,i are replaced by ^JJTi^i in all of the formulas. We get rather 
complicated homogeneous mutual distances 

2MiM2P?2 = M2(|w| + W3) 

2/^lM2P31 = + Pl)\M + {p-ivl - p.l)w3 + 2v2^pip,2Wi 

2A*1M2P23 = + Pl)\'w\ + {p.2^'1 - P-l)w3 - 2viy/ pip,2Wi. 
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In the equal mass case with mt — 1 and \w\ — 1, however, we get 

2 III 2 III 2 I I 1 

Pl2 = H + W3 P31 = \w\ + — U>i - -W3 P23 = - ^^"1 - 2"'^' 

On the other hand the Hamiltonian is 
where the norms are Euchdean. 

Example 6 (Equilateral coordinates on S^). If we use the basis of example |4] 

ei = (l,a;,w) 62 = -ei = (-1, -w, -w) 
we get simple mutual distances 

P12 = \M - Wi = \w\ + ^Wi - pis = \w\ + ^Wi + ^W2. 

CoUinear shapes form the equator W3 = with the binary collisions placed at the 
roots of unity. 

On the other hand we have a formidable conformal factor 

3mim2m3m{wl + + w^) 

In the equal mass case (m^ = 1 ) we see k — I. 

5.6. Visualizing the Shape Sphere. Having reduced the planar three-body prob- 
lem by using size and shape coordinates, we will pause to have a closer look at the 
shape sphere and the shape potential. 

Using the spherical variables w = {wi, ^2,^3) we can visualize the shape sphere 
as the round unit sphere in M'^. The equilateral basis of example [g] puts the binary 
collisions at the third roots of unity on the equator and the Lagrange equilateral 
configurations at the poles. Figure [T]) shows some of the level curves of V for two 
choices of the masses. In addition to the binary collisions shapes where — > 00 , 
there are three saddle points at the Eulerian central configurations. The Lagrange 
points are always minima of V. 




Figure 1 . Contour plot of the shape potential on the unit sphere 
wf+wl+wl — 1 in the equal mass case (left) and for masses 
TOi — l,m2 = 2,7713 — 10 (right). 
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If we use stereographic projection to map the sphere to the complex plane, we 
get the affine coordinate representation of example |4] Figure [2] shows affine contour 
plots for the same two choices of the masses. Now the collinear shapes are on the 
real axis. 




-1.5 -1.0 -0.5 0.0 05 1.0 1.5 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



Figure 2. Contour plot of the shape potential on the complex 
plane in the equal mass case (left) and for masses mi = 1, m2 = 
2,7713 = 10 (right). These plots can be viewed as a stereographic 
projections of those in figure [l] 



6. Levi-Civita Regularization 

In this section, we describe a way to simultaneously regularize all 3 binary colli- 
sion using 3 separate Levi-Civita transformations. This approach to simultaneous 
regularization was introduced by Heggie [8]. There are two versions depending on 
whether the variables Qij or the homogeneous variables Xij are used. The former 
approach was used by Heggie; we will take the latter. We begin with a review of 
Levi-Civita regularization for the Kepler problem. 

Levi-Civita showed how to regularize the two-body problem, which is to say, the 
Kepler problem. Let q d C denote the position of a planet going around an infinitely 
massive sun placed at the origin. After a normalization, the Kepler Hamiltonian is 
l|pp — Levi-Civita's transformation is the map 

z !—> = g 

together with the induced map on momenta: 

1_ 

and the time rescaling 

d _ d 
dr dt 

To understand the map on momenta, make the substitution q = va the expression 
{p, dz) for the canonical one-form. We have (p, dq) = {p, 2zdz) — {2zp, dz) which 
shows that if = 2zp then (77, dz) — {p, dq). This computation shows that the map 
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{r],z) — > {p,q) with p = j^rj, q = is a 2:1 canonical transformation away from 
the origin. Observe that r = \z\'^. Thus in terms of the new variables 

Time rescaling is equivalent to rescaling the Hamiltonian vector field. This 
rescaling can be implemented using the following "Poincare trick" . If Xh is the 
Hamiltonian vector field for H, and if ft, is a value of H, then fXn is the Hamiltonian 
vector field for the Hamiltonian H = f[H — h) provided we restrict ourselves to the 
level set {H = h}. We take f = r = \z\^ and compute that 

H=li\v\'-h\z\'-a). 

which is the Hamiltonian for a harmonic oscillator when h < 0. 

6.1. Simultaneous Regularization. Let {r,X) denote either the homogeneous- 
spherical or radial-homogeneous coordinates. To simultaneously regularize all three 

double; (;ollisions wc perform a Lcvi-Civita transformation on each of the homoge- 
neous complex variables X^j. Thus, we introduce three new complex variables 

j^ij '^ji set 

Y r2 

Define a regularizing map / : Cq ^ Cq by 

X = f{zi2,Z31,Z23) = (^12) ^31' ^23)- 

The preimage of the subspace W is the quadratic cone C with 

W: + 4i + 4a = 

and wc have f : Cq ^ Wo- Note that every X E Wo has 8 prcimages under /, 
except for the three binary collision points {Xij = some ij) which each have 4 
preimages. (Since X ^ 0, at most one of the Xij or Zij can vanish at a time on Wo 
or Co.) 

Since / is homogeneous, it induces maps fsph : — > and fpr : CP^ CP^. 
In this case we also view Zij as homogenous spherical or projective coordinates. 
These restrict to regularizing maps f^ph : S(C) S(<S) and fpr : P(C) P(W) 
where, as above, S(.),P(.) denote quotient spaces under real and complex scaling, 
respectively. 

The mutual distances become 



(56) pij = \Xij\ = \zijf 
and the mass norm is 

(57) \X(Z)\'^ = \f(z)\^ = "^l"^2Pl2 + Wimspji + TO2TO3Pi3 

mi+m2 + ms 

We will use the standard Hermitian inner product, denoted ((., .)), on 2;-space so 

(58) ||2;||^ = \z12f + 1-231 1^ + |-223|^ = P12 + P31 +P23- 

Let rjij be the conjugate momenta to Zij and let Yij the homogenous momenta 
conjugate to Xij. We extend / to a map {r,Pr,X,Y) = F{r,Pr,z,r]) by setting 

1 



PLANAR THREE-BODY PROBLEM 



33 



Then F restricts to maps T*E.+ x T*p^C^ T*M+ x T*j,^C^ and T*R+ x T*^^ 
r*M+ x T*^C^ where in (z, 7y)-variables we have the constraints re(?7, z) = for 
the sphere and {rj, z) = for the projective plane. We continue to denote these 
restricted maps by the letter F. 

The action of c £ C by translation of the momenta Yij to Yij + c pulls-back under 
F to translation of rjij by 2czij, that is, to the action 

c- (r,pr,z,'i]) = {r,pr,z,ri + 2cz). 

The momentum map for this pulled back action is 7 = + ^31 + ^23- course 
we will be interested in the level set 7 = 0. We will call this the z-translation 
symmetry of ry. 



6.1.1. Geometry ofC and the Regularized Shape Sphere. It is interesting to investi- 
gate the algebraic surface C in more detail. If we write the complex vector z Q 
BuS z = a + ib where a = re z, 6 = im z e then 



+ z: 



31 



+ ^23 



if and only if |a| 



\b\' 



0. 



This means a, b are real, orthogonal vectors of equal length = |q 
If we define a third vector c = a x & we get an orthogonal frame in ] 



= |6p = |z|V2. 
and the matrix 



(59) 



A{z) = 



ai2 ai2 C12/S 
031 hi C31/S 

023 &23 C23/S 



e 5*0(3). 



The mapping A{z) induces a diffeomorphism from the quotient space of S(C) of 
Co under positive, real scaling to 5'0(3) and hence, as is well-known, to the real 
projective space RP(3) (and to the unit tangent bundle to S"^). 

The projective curve P(C) turns out to be diffeomorphic to the two-sphere 
and, accordingly, we will call it the regularized shape sphere. One way to see this is 
to note that P(C) ~ S(C)/S^ is the quotient of S(C) under rotations. It is easy to 
see that action the rotation group on z rotates the vectors a, 6 G M'^ above in their 
own plane and leaves c = a x b invariant. It follows that the map z i— >■ c/\c\ induces 
a diffeomorphism P(C) ~ S^. 

In the sections below, we will apply the regularizing map to obtain several regu- 
larized Hamiltonians for the three-body problem. Starting with homogenous spher- 
ical variables leads to a regularized system not reduced by rotations while the radial- 
homogenous variables lead to a Hamiltonian system which is both regularized and 
reduced. In addition we will consider several way to parametrize the cone C to ob- 
tain lower-dimensional systems. Theorem [l] can be applied to show the equivalence 
of the Hamiltonian systems below, but we will omit most of the details. 



6.2. Spherical Regularization. First we will find the regularized Hamiltonian 
in spherical-homogeneous coordinates. This gives a regularization of binary colli- 
sions without reducing by the rotational symmetry. Let (r, X) be the spherical- 
homogeneous coordinates of section |4] The spherical Hamiltonian is 

H,p^,ir,pr,X,Y) = IpI + ^X(y) - -ViX). 

2 r 
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Using the formula anafogous to the one in ([7| for K{Y) and applying the regular- 
izing map gives 



nsph[r,pr,z,i])- p^-f- I + + 



(60) 



2 r2 \8mipi2P3i Srn2Pi2P23 87713^31^23 

1 / mim2 m^mi TO2"^3 



r \ P12 P31 P23 



and where 

(61) TTl = rii2Z3i - T]3iZi2 7r2 = ?7232;T2 - ?7l2^23 TTa = ?73iZ23 - r?23^31- 

Next we rescale time using the Poincare trick. One choice of time-rescaling 
factor is 1^12-231 -2^231^ = Pi2P3iP23- But since X, z are homogeneous coordinates, the 
degree-zero homogeneous function 

/„„N _ P12P31P23 _ P12P31P23 

^ ~ ipi2 + P31 + P23)^ ~ \\z\\^ 

is more appropriate. Note that by the arithmetic-geometric mean inequality we 
have < r < ^. 

The rescaled solution with energy iJ^p/i — h become the zero-energy solutions of 
Hsphir,Pr,z,ri) = T{Hsph - h). 

fr rpl |X(z)p / K1IV23 , \^2?P31 , k3ppi2 \ . , 

(63) H,^,^— + -^^-^— + ^—^^—y-W{z)-hr 

where the regularized shape potenial W is 



\Xiz)\ 

\4 



(64) W{z) = ,|_ {mim2P3iP23 + niim3pi2P23 + TO2TO3P12P31) 



Note that since z is a homogeneous variable representing [z] € S^, we have z 0. 
For a homogeneous coordinate representing a binary collision we will have exactly 
one of the variables = and ||z|| > 0. Thus H is nonsingular at these points 
and the binary collisions are regularized. 

Theorem 9. The Hamiltonian flow of H^ph on T*M.'^ x T*Cq has an invariant 
submanifold T*M.~^ x T*^^ ^^Cq defined by reirj, z) — Q and + zf^ -|- ZI3 = 0. The 
quotient of the restricted flow by scaling and by translation of rj by z represents the 
zero total momentum three-body problem with regularized binary collisions, reduced 
by translations (but not by rotations). 

The quotient space of r*p^ ^Cq by these symmetries can be identified with 
T*S(C) ~ r*RP(3). The regularizing map induces an 8-to-l branched covering 
map fsph '■ S(W) — > S(yV), that is, an 8-to-l branched covering RP3 ^ S3. The 
map is a diffeomorphism except where (exactly) one of the Zij — and Xij = 0. 
To describe the branching behavior, note that in the two-dimensional complex sub- 
space W, the set where X12 = is a complex line which corresponds to a circle 

in the sphere S(yV). The preimage of this circle will be 2 circles in the projec- 
tive space S(C). Altogether, the map is branched over 3 circles, each circle having 
pre-image 2 circles in the projective space RP'^. 
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6.2.1. Quadratic Parametrization of C . Instead of writing Hamilton's equations for 
Hsph we will describe a parametrization of the cone C which leads to a lower- 
dimensional system of equations. There is nice 2-to-l parametrization by quadratic 
polynomials which is related to the double covers of RP^ by S3, of 50(3) by the 
unit quaterions, and of SO{Z) by SU{2). 

Define a 2-to-l mapping g : ^ C C C'^ by 

(65) g : zi2^2ixiX2 Z31 ^ xj + xl Z2z=i{x\-x\) 

where Xi,X2 € C. This can be seen as a variant of a map used by Waldvogel j3D] in 
his regularization of the planar problem. But here we are applying the idea to the 
homogeneous variables X which makes it easier to blow-up triple collision later on. 

By homogeneity there is an induced map gspu ■ S'^ — ^ S{C). The induced map is 
given by the same formula except that x, z now denote homogenous coordinates for 
the points of S'^, S^. (This double covering map gives another way to see that S{C) 
is diffeomorphic to the real projective space RP'^.) The map gsph can be motivated 
in several ways. First, after omitting the factors of i, it resembles the formulas for 
parametrizing Pythagorean triples. Next, write xi = ui — i U2, X2 = M3 + i U4 and 
define the unit quaternion u = Ui+iu2+ j u^ + ku^^. Then the familiar conjugation 
map V 1—^ uvu where v is an imaginary quaternion defines a rotation R{x) on the 
3-dimensional space of w's. Up to a permutation of the columns, R{x) ~ A{z), the 



matrix of ( 59 ), and hence the conjugation map defines a map x ^-^ z. As a variation 



on this construction, define the unitary x-dependent matrix 



U ■ 



Xi X2 
-X2 Xi 



e SU{2). 



Then the adjoint representation v ^ U{x)vU(x)^^ on su{2) ~ M.'^ produces the 
same rotation R{x). 

The composition / o g^pf^ of the regularizing map and the quadratic paremetriza- 
tion gives a 16-to-l branched cover S'^ H> S'^, which becomes 8-to-l over the binary 
collisions. Each binary collision is represented by a circle in the range which has 
2 preimage circles for a total of 6 branching circles in the domain. Using stereo- 
graphic projection, it is possible to get some idea of the behavior of this remarkable, 
regularizing map. Figure |3] shows the projection of the three-sphere. The three 
transparent surfaces are tori representing the coUinear configurations with a given 
ordering of the bodies along the line. These intersect in 6 circles representing the 
binary collisions. The figure shows thin tubes around each of these circles. 

To extend 5 to a partially symplectic map G : T* M.+ xT*C'^ T*M+ xCxC^* 
we transform the momenta 7/,y so that y — rjDf{z)ov 

[yi 2/2] = [?7i2 Jysi ^723 



-2ix2 —2ixi 
2x\ 2x2 
-2ixi 2ix2 



The value of rj is not uniquely determined but any two solutions will yield equivalent 
covectors and the same transformed Hamiltonian. For example, we could take 

'712 = 

G restricts to G : T*R+ x T*p,^C^ T*R+ x T*^^,^^'^^* where T*p,,€? = {{x,y) : 
Te{y,x) = 0} and T^^^^cC'* - {(z^v) ■ z € C,re(??,4 = 0}. 









^ ( yi 


_ yA 


'731 = I 






??23 = T — 








X2J 


4 \xi 


X2J 
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The regularized spherical Hamiltonian becomes 
(66) 

^-ph- 2 + r2|la:||in 256mi 256m2 ^ 2567713 ; r^ ' 
T^'i ^ yiX2 + y2Xi ^ yiX2 - y2Xi tt's = yixi - 1/2X2 

10 |2 I 2 I 2|2 I 2 2|2 

II ||2 oil l|4 ,1 lYf m2 mim2pf2 + ™i™3Pii + m2m3p^3 

IfII =2||a;|| =pi2 + P3i+P23 \X{x)\ = ■ ■ . 

mi + TO2 + TTT-a 

Note that H is invariant under the scaling symmetry {x,y) — )■ {kx,k~^y), k > 0. 
The corresponding Hamiltonian system on the 10-dimensional space T*(IR"'' x C^) 
can be reduced to the expected 8 dimensions by restricting to the invariant set 
r*]R+ X T*pifi^ and then passing to the quotient space under scaling. 

6.3. Projective Regularization. Next we will get a regularized version of the 
reduced three-body problem. Let (r, X) be the radial-homogeneous coordinates of 
section [5] For a fixed angular momentum, we have the reduced Hamiltonian on 

H,{r,p^,X,Z) = + ^) + ^-^KiZ) lv{[X]). 
After making the Levi-Civita transformations, fixing an energy and changing timescale 



by the factor r from ( 62 1 we obtain a regularized reduced Hamiltonian 



rp2 rfi^ , |X(z)|2 f\ni\^p23 , k2pP3i , k3pPi2\ 1, 

r 



(67) H, = — + ^ + -^^-^ + ^— + ^—]~-W{i)-hr 



where the various quantities appearing in the formula are given by (56), (57), 



(58) and (64). The only difference between the spherical and projective Hamiltoni- 
ans is the term involving /x^. We also impose the extra constraint im(?7, z) — Q and 
there will be extra curvature terms in the differential equations. 

To find the curvature terms we need to pull-back the Fubini-Study form under 
the regularizing map X = f{z),Xij = zf,. The Fubini-Study metric on z-space is 



derived from the standard Hermitian metric on C by a formula analogous to (31 ). 
We can express its restriction to C in terms of a tangent vector field as we did in 
lemma [2] The analogous formula to ( 32 ) is 



(68) {{V, W))FSfi = V,W <^ TxS 

where e(z) is a Fubini-Study unit vectorfield tangent to C and normal to z. For 
example, observe that if z S Cq = C \ then the vectors z, z, T form a Hermitian- 
orthogonal complex basis for TJZ'^ where 

(69) T = Z X Z = (Z31Z23 - Z23Z3i,Z23Zi2 - Zi2Z23,Zi2Z3i - Z31Z12). 

Hence we can take e = p||T = (z x z)/||z||. This gives 

„ „ S(g)S 



(70) 



IFS,C 



WAV 
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where S is given by any of the following formulas 

^ _ ((z X z,dz)) _ ||z||(zi2(iz3i - 2:31^212) 

(71) 



\\'A\ Z2Z 
Z\\{z23dzi2 - 2:12^^23) ||2||(23lC?223 - 2231^231) 



231 



212 



For example, the first version is just S = ((e, dz)) and the second is obtained by 
eliminating Z23, dz23 using the equations 2I3 = —zf2~ z'^i and 223^223 — —212^212 — 
231 '^231. 

Using these formulas, we find that the pull-back of the Fubini-Study metric on 
W is a conformal multiple of the Fubini-Study metric on C 

Lemma 4. The pull-back of the Fubini-Study metric on W is given by 

/*(•, ■)fs,w = A(2)((., .))fs,c 
where the conformal factor is 

)^ _ 4toiTO2"^3 P12P311O23II2IP _ 4mmim2m3{pi2 + psi + P23) P12/O31P23 
m|X(z)|4 (miTO2/0^2 + fni'nT'SPli + m2m3pl^y 

and where pij — \zij\'^ . 



Proof. Equation (34) shows that we need to compute the puUback f*a, where a is 
given by ( 35 1 . Using the first formula for <t gives 



f*a = 22i223l(i23i - 223i2i2(i231 = 2212231223^ 



Hence 



TOiTO2™3 



4toiTO2™3 



rn|A(2jp m|A(2)p 
Now use (571, ( 58 ) and ( 70 ) to get the formula in the proposition. QED 

Similarly we can pull-back the Fubini-Study cometric on W and compare it with 
the dual Fubini-Study metric on C. The formula analogous to (33) is 

(73) 



l'7llkc = l(^,e)| 



2 _ \{V,ZX Z)|2 



l^ll 



z,pr 



This is a degenerate quadratic form, invariant under z-translation of 77, which rep- 
resents the Fubini-Study cometric on C. 

The next lemma relates this to the pull-back of the Fubini-Study cometric on W 
and hence, to the shape kinetic energy. 

Lemma 5. The pull-back of the Fubini-Study cometric on W is 



F* 



\FS,W 



= X- 



|2 

\fs,c 



where A is given by (72). Hence the shape kinetic energy in regularized coordinates 
is 

1,-1|,||2 _ 1 liV^ZX 2)P 
2^ ™PS.C-2 ;^||^||2 • 
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Proof. Equation (36) shows that we need to compute the puhback F*a, where a 
is given by (37). Using the second formula for a gives 

|^23p^*^ _ (??31^12 - ?7l2^3l)^:23 
I^P 2zi2Z3iZ23 

and there are two similar equations from the third and fourth formulas. Adding 
these gives 

Therefore, 

^*||,|2 m\X{z)n{r,,zxz)\' ^ m\X{z)\^ , 

-'^■^'^ 4mim2r7l3Pl2P3lP23||2||'^ 4mlTO2"^3Pl2/^3lP23||^||^ 

Comparing with the formula for A completes the proof. QED 

It follows from the lemma that we have an equivalent reduced, regularized Hamil- 
tonian 

^ rp^ m|X(z)|4|(?7,z x z)p _ 1 _ 

2 2r^ 8mim2m3r2||z||io r 

The factor of A in the Fubini-Study two-form and the factor of A^^ in the shape 
kinetic energy cancel out in the interior product defining the curvature term. Re- 
membering the timescale factor r we find that the the curvature term is 

(75) Tcurv — i V 

which is added to the right hand side (i.e. to — ^) of the Hamilton's equation for 
ri- 

Theorem 10. The Hamiltonian flow of on r*M+ x T*Cq has an invariant set 
r*M+ X T*^(,C^ where {rj, z) = and -|- z^^ + = with symplectic structure 
given by the restriction of the standard form minus 2iiXflps. The quotient of the 
restricted flow by the complex scaling symmetry and by z -translations of rj represents 
the three-body problem with zero total momentum and angular momentum /i, with 
regularized binary collisions, reduced by translations and rotations. 

The regularized, reduced Hamiltonian iJ^, together with the curvature term 
gives a system of differential equations on the 14-dimensional space T*(M"'" x C^) 
with variables {r,pr, z,ri). The six-dimensional quotient space of r*M+ x T*^.qC^ 
is diffeomorphic to T*M+ x T*P(C). Instead of writing these 14-dimensional differ- 
ential equations we will describe several ways to parametrize the regularized shape 
sphere P{C) to arrive at lower-dimensional systems of equations. 

6.3.1. Quadratic Parametrization of the Regularized Shape Sphere. We can param- 
etrize C using the same quadratic map g : — !■ C C C'^ 

Zi2 = 2ixiX2 Z31 = x\ + x\ Z23 = - xl) 

as in section |6.2.1[ Since g is homogeneous with respect to complex scaling, it 
induces a map gpr : CP^ — >■ P{C) from the projective line onto P{C). Although g 



and the induced map gsph of S in section 6.2.1 are both 2-to-l, the extra quotienting 
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makes gpr a diffeomorphism. This shows again that P{C) is diffconiorphic to the 
two-sphere. The same partiahy symplectic extension G : T*E+ x T*C^ T*R+ x 
C X C^* restricts to a map G : T*R+ x T;^C^ T*R+ x T*^ c'^^ where T^^C^ = 
{ix,y) : {y,x) = 0} and T;^^^'^^ = {(2,77) : z G C, (77,2) = 0}. ' 

If we use ( 74 ) together with the formula ([73| for the dual Fubini-Study metric 



we obtain, after some simplification, the reduced, regularized Hamiltonian 
= ^ + ^ + ^lyi^2 - x,y2\' - Iwix) - hr 



(76) 



W{x) = ^f^u)} ("117712^31^23 + rnim3Pi2P23 + TO2"l3Pl2P3l!l 



P12 = |2a;ia;2p p3i = |xi+X2p P23 



on ||4 , , lYf M2 TOiTO2P?2 + "^i"^3Pii + nT-2"l3Pi3 

2||a:;|| =/5i2+P3i+P23 = 



7711 + 7712 + ni3 



We have 
(77) 



P12P31P23 T _ 77i|X(a;)|'' 



8||x||-'-^ A 64toi?7721'13 "'"^ 

We have the complex constraint (y, a:) = and the system is invariant under 
complex scaling symmetry {x,y) — > {kx,y/k), k € Co- Applying the constraint 
and passing to the quotient space reduces the dimension from 10 to 6. As usual, 
Hamilton's differential equations will have a curvature term 



j,2 *y 



added to the y equation. 



6.3.2. Dynamics in regularized affine coordinates. As in section [5^ we can use afhne 
local coordinates on CP^. Every projective point [xi,X2] £ CP with xi ^ has a 
representative of the form 

[a;i,X2] = [l,z] = [l,x + iy] 

where x,y ^M.. The appropriate momentum substitution is 

yi = -zC y2 = C- 

where C. = a + i/3 Cz C* is a momentum vector dual to z. 

We get a Hamiltonian system with 6 degrees of freedom 
(78) 



W{x, y) = ^'^^^'^^\sf, ("ll"l2P3lP23 + "ll"l3Pl2P23 + "l2"l3Pl2P3l) 
(1 + + 7/^)° 

P12 = 4(a;2 + 7/2) = 7/2)2+4^2^2 ^ (l_a;2+ ^2)2+4^2^2 

\\xr^l + x'+ 7/2 \X{x, 7/)|2 = ^^i"^2P^2 + m^m^pj, + 77127713/723 

77li + 7II2 + 7II3 

The Fubini-Study form becomes 

dx A dy 



Qps — 



(1 + a;2 7/2)2 
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and 



P12P31P23 



8(1 + a;2 + y2-)6 X 64mim2m3 (1 + a;2 
Hamilton's equations with the curvature term are 
r 



Pr 



(79) 



rpr 

T (1 + : 
A 

T (1 + : 
A 



(i + iznici^+v 



1 



W^(x,y) 



2^2 



y^f 



^ ,,2l 



r 



/3= ^[-Ay-2T^a-^/i^] 



9 

pi 



-h] 
h] 



where 



2A 



6.3.3. Dynamics in regularized spherical coordinates. Instead of using projective or 
local afSne coordinates, one can map the regularized shape sphere to the unit sphere 
in a particularly elegant way to do this is to use the diffeomorphism between 



A{z) 



C and 50(3) described in section 6.1.1 

Given z e C we write z ^ a + ih where a, 6 € M'^ and then define c = a x 6 e M''. 
We saw that the matrix 

ai2 ai2 C12/S 
: 031 ^31 C31/S 
[023 ^23 C23/s_ 

is in 50(3), where = \z\^ /2 = |a|2 = \b\^ = \c\ 

We will work homogeneously and define a map 5 : C — J' M^^ 

g{z) — c — re(z) x im(z). 

By homogeneity, there is an induced map gpr '■ P(C) — > S(M'^) ~ S'^ where we view z 
and c as homogeneous coordinates with respect to complex and positive real scaling 
respectively. 

The orthogonality of the matrix A{z) can be used to derive some useful formulas. 
Since the rows as well as the columns are unit vectors we find 



Pij =a-j +h% = 



^23 



„2 
-12 



which gives the beautiful formulas 

/or\\ „ ^31 + ^23 „ '^12 + ^23 „ CJ2 + eg 

(o'J) P12 = — n — P31 = — n — ^23 = — n — 

\c\ \c\ \c\ 

for the homogeneous mutual distances. Similar formulas were given by Lemaitre 



^-31 



Next, consider the quantity 
Z12Z31 = 012031 + 612631 + i (012631 - 031612) = (ai2, 612) • (031, 631) + i C23. 
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Using the orthogonality of the rows we can express this entirely in terms of c. We 
find 

C12C31 . _ C23C12 , . _ C31C23 , . 

-2122:31 = n ^ * "^23 2:23-212 = n 1" * ^31 -231^23 = n 1- « C12. 

\c\ \c\ |c| 

These last formulas allow us to write down local inverses for Qpr- Namely consider 
the map /112 : ^ 

/112(c) = \c\zi2{zi2,Z3i,Z23) = |c|(2l22l2, 212-231, ^12^23) 

= (C31 + cls, -C12C31 + i \c\c23, -C12C23 - i \c\c31). 

If 2^12 7^ 0, then /112(c) represents the same projective point in P(C) as z does so 
/112(c) give a local inverse for the projective map gp^. There are similar partial 
inverses /131, /i23- 

To find the regularized, reduced Hamiltonian system we need to convert the 
Fubini-Study metric and its dual norm (i.e. cometric) to c-coordinates. The spher- 
ical analogue of the Fubini-Study metric is the spherical metric 

\c\^{dc,dc) — {dc,c){c,dc) |c x rfc|^ 

-/^P^ = ^4 = |c|4 

where we are using the Euclidean inner product on R^. We will see that 

*/ \ o// \\ 2|((z X z, dz))\'^ 
9 {; -Isph = ^{{■, -nFSfi = ii^^jje • 

To see this, note that z x z = —2i ax b = —2i c. Hence 

1 

dc = -{dz X z + z X dz) . 

This, together with the fact that {{z, z)) = on C leads, after some algebra, to the 
pull-back formula. Correspondingly, the Euclidean solid angle form pulls back to 
twice the Fubini-Study form, hence 

Ar^FS.c = 9* ttA? (cic;c2 a dc3 + C2dc3 A dc\ + C3dci A dc2). 
2|cr 

Let 7 S R^* be a dual momentum vector to c G M'^. From the spherical scaling 
we will have 7 • c = 0. If we split the momentum vector r] into real and imaginary 
parts, 7] = u + iv then the momenta transform via 

U ' C V ' c 

|cr |cr 

Prom this we find that the dual spherical norm \^\lpf^ = I7 x cp = |cp|7p corre- 
sponds to lll-lll'sc- reduced, regularized Hamiltonian 



\X{c)\ 
(81) ■ 8|c 



^i^) = ol„l6 ("^l"^2P3lP23 + mim3Pi2P23 + W2m3Pi2P3l) 



P12 = C^i + CI3 P31 = CI2 + CI3 P23 = C?2 + C31 

oui2 ^ , . , ^ \y(^W^ rnim2pl2 + mmispli + m2W3pi3 

^C =Pl2+P31+P23 = ; ; 

mi +1712+ 1713 
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We have 

/O12P31P23 T m\X{c)\'^ 



%\c\^ A 64mim2m3 |c|8 ■ 

Here we have redefined pij to ehrainate the factors of |c| and placed these factors 
elsewhere in the formulas. The curvature term is 

, - 2^r 

(82) TcuTv = rr^ 7 X c- 

6.4. Visualizing the Regularized Shape Sphere - LeMaitre's Conformal 

Map. The map of projective curves fpr '■ P(C) P(yV), induced by the squaring 
map, can be visualized as a map of the two-sphere into itself. Indeed this is the 
point of view taken by Lemaitre in [13 . 

The map is a four-to-one branched covering map with octahedral symmetry (see 
figure |4|. The map is generically four-to-one except at the binary collision points 
where it is two-to-one. In the figure, each octant of the regularized sphere maps 
to one or the other hemisphere of the unregularized sphere. Thus, for example, 
the north pole of the unregularized sphere (representing a Lagrangian, equilateral 
central configuration) has four preimages which lie in alternate octants. Each binary 
collision point on the equator of the unregularized shape sphere, has two preimages, 
which lie on a coordinate axes of the regularized sphere. 

The three-dimensional sphere of figure [3] is just the preimage of the regularized 
two-sphere sphere in figure |4] under a Hopf-map. Each point of the two-sphere 
determines a circle in the three-sphere. The three large tori in figure |3] are the 
preimages of the collinear circles in the two-sphere (where the coordinate planes 
cut the sphere). The six tubes in figure [s] are the preimages of small circles around 
the binary collision points (where the coordinate axes cut the sphere). 

7. Blowing Up Triple Collision 

Our systematic use of the radial coordinate r together with the homogeneous co- 
ordinates used to describe the shape make it easy to implement McGehee's method 
for blowing-up total collision. We need only rescale momenta and change the 
timescale. The changes can be made before or after reduction. The changes are 
non-canonical so destroy the Hamiltonian character of the equations. We will de- 
scribe the general method for the rotation-reduced and unreduced cases and then 
make some comments on the results of applying it to some of the Hamiltonians 
described above. 

7.1. Before Reduction. Consider a Hamiltonian of the general form 

(83) H{r,Pr,X,Y) = ^B(X){Y,Y) - -V{X) + [\a{X)pI ~ C(X)] 

2r^ r 2 

when expanded in powers of r. This covers the rotation-unreduced Hamiltonian 
Hsph of section [4] and the corresponding regularized Hamiltonians Hsph{i~,Pr, z,r]) 



and Hsph{r,pr,x,y) of section 6.2 (after changing the names of the variables). For 
the unregularized Hamiltonian H^ph we have 

A(X) = 1 C{X) = 

while for the regularized Hamiltonians Hsph we have 

A{X):^t{X) C{X)^hT{X). 
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The quantity B{X){Y,Y) represents the non-radial part of the kinetic energy. It 
is a quadratic form in Y which we represent by a symmetric matrix B{X) depend- 
ing on X. The dependence of i? on X must also be quadratic since H must be 
homogeneous of degree with respect to the scaling {X, Y) i-)- {kX, -^Y). 

Let /(r) be a positive, real-valued function. We will introduce a new timescale 
such that ' — fir)' ■ The usual choice is McGehee's scaling factor /i(r) — but 

we will also consider /2(r) = (^^qif) ^ which has better behavior for large r. (With 

the 1st choice solutions can reach r = oo in finite time.) For any such /(r) we 
replace {pr,Y) by rescaled momentum variables 

(84) 

The shape variable X remains the same. When we make these substitutions of inde- 
pendent and dependent variables in the Hamilton's differential equations resulting 
from ( 83 1 we get 
(85) 

r' = A{X)vr 

v' = + r{\nv)r)A{X)v^ + B{X){a,a) - v{r)V{X) 
X' ^ B{X)a 



r{\nv)r)A{X)va ~ \ Ax - \Bx{a,a) + v{r)Vx + rv[r)Cx 



2' V V ' 2 2 

where ;/(r) = f{r)'^/T^ and where the subscripts denote differentiation. For McGe- 
hee's scaling /(r) = /i(r) = we have i'(r) — \,i\ni')r = and the equations 
simplify considerably. For /2(r) we have v{r) = (1 + r)"'^ and both v and (Inz^)^ 
are still smooth all the way down to r = 0. 

Writing the energy equations Hgph — h or Hsph = in terms of the rescaled 
momenta gives 

(86) \a{X)v^ + \B{X){a, a) - v{r)V{X) = ru{r)C{X). 

For example if use the rs rescaling with Hsph, we have 

A=l B{X)^\X\'^Bq C = V{X)^\X\^ 



where Bq is the constant symmetric matrix We get the blown- up differential 
equations 

r' = vr 

v' ^^^v-" -\X\''Bo{a,a) + V{X) 
X' = \X\^BQa 

a' — ^2^*^ ^ -So(a, a)X + Vx 



with the energy relation 



+ ]^Bo{X){a,a) ^V{X) ^ rh. 
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The regularized equations arising from Hgph are considerably more complicated 
due to the B{X) terms (or rather the B{z) or B{x) terms.) Instead of writing them 
explicitly we will just ma ke some observations about them. Consider, for example, 



Hsph{f,Pr,x,y) from (66). B(x) will be a complicated, 4 x 4 real matrix arising from 
the second term in ( 66 ). The phase space before blow-up is r*M+ x T*C^ ~ (0, oo) x 
K X X C^. In addition to the energy relation Hgph = 0, we have re(y, x) = and 
the scaling symmetry by positive real numbers so there is an induced flow on an 
quotient manifold of real dimension 7. After blow-up we have variables {r,v,x,a) S 
[0,oo) X M X X C^, where we have extended the flow to the collision manifold 
where r — 0, which is an invariant set for the differential equations. We have a 
real-analytic vectorfield on this manifold-with-boundary. Imposing the constraints 
and passing to the quotient under scaling gives a real-analytic vectorfield on an 7- 
dimensional manifold-with-boundary representing the planar three-body problem 
on a fixed energy manifold, with all binary collisions regularized and with triple 
collision blown-up. Note, in particular that the regularization of binary collisions 
passes smoothly to the boundary. 

We claim that if the timescale factor /(r) = f2{r) = (r/(r + l)) 3 is used, then the 
differential equations define a complete flow on [0, oo) x M x x and hence the 
induced 7-dimensional flow is complete. Since the differential equations are smooth, 
the only obstruction to completeness would be orbits which become unbounded in 
finite time. It is well-known that, with the usual timescale, such orbits do not exist 
for the three-body problem. It follows that if we use only bounded time-rescaling 
factors, the same will hold for the modified differential equations. McGehee's factor 
r2 is unbounded and it is possible for orbits to escape in finite time. Indeed, there 
are solutions of the three body problem for which r{t) = 0{t) as t oo with 
respect to the usual time-scale and these will reach infinity in finite rescaled time. 
The factor /2, while producing less elegant differential equations, eliminates this 
problem. 

7.2. After Reduction. The rotation-reduced Hamiltonians and their many 
regularized forms i?^ have the general form 

(87) H,,{r,pr, X, Z) = ^[B{X){Z, Z)+A{X)fi^] - ^V{X) + [^A(X)pl - C{X)] 

(after changing the names of the variables). The only new term here, when com- 
pared to the Hamiltonian of section |7.1[ is the quadratic term in the angular mo- 
mentum fi. We have a momentum constraint {Z, X) = and there will be a 
curvature term 



added to the Z equation. As in section 7.1 the unregularized Hamiltonians 7?^ we 
have 

A{X) = 1 C{X) = 

while for the regularized Hamiltonians we have 

A{X)^t{X) C{X)^hT{X). 

As in the last section, the variables X, Z can denote either homogeneous coordinates 
on the cotangent bundle of projective space, before or after Levi-Civita transfor- 
mation, or they can be local holomorphic coordinates on the cotangent bundle of 
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the shape sphere or of the regularized shape sphere P(C) (see the examples below). 
Our computations immediately below hold for all these cases. 

We rescale time and the momenta as in equation (84) with Z replacing Y . We 
must also rescale angular momentum according to 

h or = become 
(89) ^^A{X){v' + li^) + ^B(X)(a,a) - v{r)V{X) = rv{r)C{X) 

where 
(90) 

so that = 1 for / = r3 and i/ = (1 + r)^'^ for / = /2 

In order to express the differential equations succinctly, let 



Then energy equations iJ^ 
1 



-3 



1 



i^=-A(X)(«2+^2) 



1 



+ -S(X)(a,a) 



denote the blown-up kinetic energy and let 

(91) 0(r) = 

Then the equations of motion are: 
" = A{X)vr 
= <i){r)A{X)v'^ + 2k- v{r)V 
= (l){r)A{X)vij, 
= B{X)a 
= (t){r)A{X)va 



(92) 



Kx + v{r)Vx+ru{r)Cx+T,, 



where 



= —2ifla or — 2iflT{X)a 
for the unregularized and regularized cases, respectively, 
equation can also be written 

v' = i(j)+ l)AiX)v^ + B{X){a, a) + A{X)jl^ ■ 

In these equations, we are regarding /i as a new variable subject, by definition, to 
the constraint 



We remark that the v' 



v{r)V{X). 



(93) yfM-TKOM 

where /i is the old angular momentum constant. This point of view is necessary to 
make the curvature term smooth at r = 0. 

As in section [7rH all functions of r extend smoothly to r = 0. If we start with one 
of the regularized Hamiltonians then for the resulting differential equations, all 
binary collisions have been regularized and the triple collision blown-up. We obtain 
a fiow on a manifold-with-boundary of dimension 5 after fixing /i, setting iJ^ — 0, 
imposing the constraint on /i, the constraints that X ^ C and (Z, X) — {) and 
passing to the quotient under complex scaling. Binary collisions are regularized 
for all values of ji and if the time rescaling is done using f2{'r), the fiows on these 
manifolds will be complete. 
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It is well-known that triple collisions are possible in the three-body problem only 



when /i = 0. In this case, equation (93) shows that either p, — or r — 0. Both of 
these submanifolds are invariant sets for the dynamical system. The 5-dimensional 
manifold-with-boundary with the above constraints and with /2 = represents 
the closure of zero-angular-momentum three-body problem. The 4-diniensional 
manifold where fl — r — forms the boundary. Even though orbit with fj, ^ 
cannot have p — > 0, the part of the collision manifold {r — 0} where /i ^ is 
relevant for studying low-angular momentum orbits passing close to triple collision 

We will now present a couple of versions of the regularized, reduced and blown-up 
differential equations for the three-body problem. 



Example 7 (The blown-up regularized affine equations). In section 6.3.2 we used 
affine local coordinates on the regularized shape sphere to obtain a regularized 
Hamiltonian H(z, C) with 6-degrees of freedom. (We wrote z~x + iy,(^ — a + if3 
in section 



6.3.2 ) Comparing with the general form (87) we have 



AiX)=riz) BiX)iZ,Z) = ^{l + \z\Y\C\' 
C{X) ^ hT{z) V{X) ^ W{z). 
For convenience, we recall from section [6. 3. 2| that 

Pl2=4|z|2 p3i = |l + z2|2 p23^|l_;22|2 

^(^) = /j^l^mfi ('7^lTO2/53lP23 + rnim3Pi2P23 + m2m3Pi2P3l) 
(1 + Izl-^jb 

I vf m2 _ mim,2pl2 + minispli + m2m3pl3 

- \ , 

nil +m2 + ma 

_ P12P31/023 _ ■m\X{z)\'^ 

~ 8(1 + |z|2)6 A ~ 64mim2m3 (1 + |z|2)8 ■ 

As per the preceding subsection, we continue to write the rescaled momentum 
variable as a (thus: a = {f /r'^)C) trusting that there will be no confusing with the 
previous use of a. The rescaled kinetic energy satisfies 

2K = Tv^+rij:' + U'^ + \z\YH^- 
A 

Then the regularized, blown-up equations read: 
r' — T{z)vr 

v' = (j){r)T{z)v^ + 2K - i^{r)W{z) 
(94) p' = (l){r)T{z)vfj. 

z' = (l + |z|2)2„ 

a' — (f){r)T(z)va — + v{r)Wz + rv{r)hTz{z) — 2ipT{z)a 
The possibilities for vir), (j){r) are described in the previous subsection, in equations 



(90 911 



We have 7 variables, (r, u, /i, z, a) £ [0, oo) xMxMxCxC. The constraints are 

1 . ^ , T It 
2 A^ 



-t{z){v' + P^) + ^^{1 + \z\^f lap - v{r)W{z) = rv{r)T{z)h 
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and 

Example 8 (The blown-up regularized spherical equations). In section 6.3.3 



we 



used spherical-homogeneous variables c = (ci, C2, C3) to give a global representation 
of the regularized shape sphere. We found a regularized Hamiltonian Hf^{r, c,pr,"f). 
Comparing with the general form (|87[) we have 



A(X) = r(c) BiX){Z,Z) = 2-\c\^\j\' 
C{X) = hT{c) V{X) = W{c). 
We recall the formulas 

P12 = C^l + C23 P31 = C?2 + C23 P23 = C12 + C31 
W{c) = ^-^rr^ {mim2P3lP23 + "lim3Pi2/923 + "l2m3Pi2P3l) 



\X{c) 



mim2p\2 + mimzpli + m2m3,pl^ 

m,i + m2 + 
P12P31P23 T _ m\X{c)\'^ 



8|c|^ A 6477117712^3 |c|^ 

With a = (//7'^)7, the rescaled kinetic energy satisfies 

2k ^Tv^ + tP? + 2Uc\^\a\^. 

A 

Then the regularized, blown-up equations read: 
r' — T{c)vr 

v' = (/)(r)r(c)7;2 + 2K- iy{r)W{c) 
jj! = (j){r)Tvji 

^^'aM'^I" 

a' = (/)(7')T(c)7;a — Kc + i^Wc + rv{r)hTc{c) H — a x c 

|c| 

We have 9 variables, (r, 7;,/i,c, a) G [0,oo) x M x M x Mq x M?. However, (c, a) 
are homogeneous variables. They satisfy (a, c) = and the equations are invariant 
under the scaling (c, a) — (fee, ija). Taking this into account, we have an induced 
system on the 7-dimensional quotient space [0,oo) x M x M x T*S^. The energy 
and angular momentum constraints are 

(96) \t{c){v^ + f?) + \\c\^ \a\^ - v{r)W{c) = rv{r)T{c)h 

2 A 

and 

Vrp = \/i4r)p 

giving a subvariety of dimension 5. 

A nice alternative to the quotient construction is just to observe that (a, c) — 
implies that |c| is invariant under the differential equations (95). Instead of 
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quotienting by the scaling symmetry we can simply restrict c to the unit sphere. 
Let 



M{h, fi) = {{r, V, fi, c, a) : \c\ = 1, {a, c) = 0, /i = 



^,(96) holds }. 



Then Ai{h, /i) is a 5-dimensional submanifold (or subvariety when /z = 0) of [0, oo) x 
E X E X Mj^ X which is invariant under (95 1. The flow on M{h,^) globally 
represents the planar three-body problem reduced by translations and rotations, 
with all binary collisions regularized and with triple collision blown-up. 
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Figure 3. Stereographic projection of showing the preimage 
under the regularizing map of the colhnear configurations and smaU 
tubes around the binary coUision circles. 
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Figure 4. The regularizing map is a four-to-one branched cover 
of the two-sphere with octahedral symmetry. Each octant of the 
regularized shape sphere (left) maps onto a hemisphere of the un- 
regularized shape sphere (right). The planes represent coUinear 
configurations. The figure also shows level curves of the unregu- 
larized shape potential and their preimages in the equal mass case. 



